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Abstract. 

We recall various multiple integrals with one parameter, related to the 
isotropic square Ising model, and corresponding, respectively, to the n-particle 
contributions of the magnetic susceptibility, to the (lattice) form factors, to 
the two-point correlation functions and to their A-extensions. The univariate 
analytic functions defined by these integrals are holonomic and even G-functions: 
they satisfy Fuchsian linear differential equations with polynomial coefficients 
and have some arithmetic properties. We recall the explicit forms, found in 
previous work, of these Fuchsian equations, as well as their russian-doU and direct 
sum structures. These differential operators are very selected Fuchsian linear 
differential operators, and their remarkable properties have a deep geometrical 
origin: they are all globally nilpotent, or, sometimes, even have zero p-curvature. 
We also display miscellaneous examples of globally nilpotent operators emerging 
from enumerative combinatorics problems for which no integral representation is 
yet known. Focusing on the factorised parts of all these operators, we find out 
that the global nilpotence of the factors (resp. p-curvature nullity) corresponds 
to a set of selected structures of algebraic geometry: elliptic curves, modular 
curves, curves of genus five, six, . . . , and even a remarkable weight-1 modular form 
emerging in the three-particle contribution x'^' of tli6 magnetic susceptibility of 
the square Ising model. Noticeably, this associated weight-1 modular form is also 
seen in the factors of the differential operator for another n-fold integral of the 
Ising class, , for the staircase polygons counting, and in Apery's study of C(3). 
G-functions naturally occur as solutions of globally nilpotent operators. In the 
case where we do not have G-functions, but Hamburger functions (one irregular 
singularity at or oo) that correspond to the confluence of singularities in the 
scaling limit, the p-curvature is also found to verify new structures associated 
with simple deformations of the nilpotent property. 
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1. Introduction 

Generating large series expansions of physical quantities that are quite often defined 
as n-fold integrals is the bread and butter of lattice statistical mechanics, enumerative 
combinatorics, and more generally theoretical physics. The n-fold integrals considered 
in theoretical physics are integrals of some more or less simple algebraic integrands: 
they are therefore holonomic pQ [5] . We actually found explicitely [31 13] the (highly 
non-trivial) Fuchsian linear ODEs satisfied by the first n-particle contribution x'"' 
of the magnetic susceptibility of the isotropic square Ising model for n = 3,4 (and 
n = 5 modulo a prime [5]). Mathematicians use to say of such n-fold integrals of 
algebraic integrands that they are derived from geometrif (DFG), which means that 
they can be interpreted as period^ of some algebraic variety [5] closely related to 
the algebraic integrancJl. Leaving behind all the cohomology that can be done on 
these algebraic varieties and other mixed Hodge structures [TOj, let us just remind 
that such DFG quantities are remarkably selected and structured. Considering 
the roots of the indicial polynomials of these ODEs (the critical exponents), one 
finds out that all the critical exponents of all the singularities of the corresponding 
minimal order Fuchsian linear ODEs are necessarily rational numbers [U lllj . Coming 
back to series expansions, these DFG n-fold integrals necessarily correspond to, 
not only Gevrey series, but convergent series and, often, to arithmetic Gevrey 
series [HI [121 HSl [HI [TTl [151 [16] . There is a notion of order of Gevrey series: order zero 
corresponding to G-functions [TTl [T71 [TSl UHl [22] , and order one to "Hamburger" [^T] 
functionE[§|, that is to say ODEs with an irregular singularity onl^ at oo. These n-fold 
integrals, corresponding to Fuchsian linear ODEs, are thus necessarily G-functions i.e. 
solutions of linear differential equations with arithmetic properties [HI [HI [19] . In a 
series of papers [Ml [25l [26l [23 [28l [291 [30] the Chudnovskys underlined the crucial 
role of this fundamental class of functionijj. They proved in [26] that solutions of 
linear differential equations satisfying an arithmetic growth property, the G-propert}|3, 
have special geometric properties in the sense that the corresponding minimal linear 
differential operators are globally nilpoten^^ From an arithmetic, as well as effective 
(i.e. computational) view-point, globally nilpotent linear differential operators [32|,l33j 

II See Picard-Fuchs operators \E\ and Gauss-Manin connection jTj. 

t These n-fold integrals can also be seen as the "diagonal" of an algebraic expression closely linked 
to the algebraic integrand [9]. 

§ In that case the integrand in these ra-fold integrals is not an algebraic function anymore. 
^ The irregular singularity is at oo or 0, but not and oo. Mathematicians do not like to consider 
direct sums of ODEs corresponding to G-functions and "Hamburger" functions 1211 : the sum of 
Gevrey series of different order is not summable and requires the multi-summability introduced by 
Ecalle l22l[23l . 

t First introduced by Siegel [31) . 

tl The G-property is an arithmetic growth property, on the coefficients of a solution-series. 

ft In particular, all solutions of such globally nilpotent operators with algebraic initial conditions are 

G-functions. 
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correspond to highly selected structures with a large number of remarkable properties. 
In particular, their wronskians are A'^-th roots of rational functions. In fact, this 
property holds as soon as the indicial polynomials for each singularity have integer 
coefficient^. A much more selected property is that all the critical exponents of all the 
regular singularities of these ODEs are necessarily rational numberJ^ Recall, however, 
that the rationality of critical exponents is a consequence of the DFG structure. 

Unfortunately this Grothendieck's geometry view-point is not very well-known in 
physics, and thus, like Monsieur JourdairQ, theoretical physicists often study series 
expansions (low or high temperature series expansions, generating functions, . . . ), or 
divergent series, without knowing that they are Gevrey series (and arithmetic Gevrey 
series) and often G-functions [l8l|T9]. In particular, they take for granted the rational 
character of the critical exponents, or the algebraic simplicity of the wronskians of the 
ODEs they encounter. Many simple remarkable results on form factors, or non-trivial 
identities on some well-poised hypergeometric series and other Bayley pairs [371 13H] , 
are not sufficiently recognized as a straight consequence of the fact that G-functions 
naturally occur. Along this hypergeometric linqjj, the paradigm of functions that 
can be interpreted as periods on an algebraic variety are the hypergeometric function 
2F1 (or sometimes complete elliptic integrals of the first or second kind K, E), or, 
more generally, hypergeometric functions n+iFn for some selected arguments. Dwork 
conjectured [ID] that globally nilpotent linear differential operators of second order, 
are necessarily reducible, up to a rational pull-back and up to the N-th root of a 
rational function, to the hypergeometric functions 2F1. This initial conjecture was 
ruled-out by KrammeJ^ who provided a counterexample |41j which comes from the 
periods of a family of abelian surfaces over a Shimura curve ( "wrong" elliptic curves) 
\ {0, 1, 81, 00}. Later other examples that are not even associated with arithmetic 
Fuchsian lattices, or Shimura curves, were also found [42]. Now the conjecture is 
rephrased to embed hypergeometric functions 2F1 and such counterexamples. In 
this paper, we will call "Dworkian" a globally nilpotent linear differential operator 
of second order corresponding to this "extended" conjecture. 

We will try to promote a DFG Grothendieck's viewpoint using a learn-by- example 
approach which focuses on the (quite arithmetic) notion of global nilpotence of the 
linear differential operators of various holonomic quantities we already encountered in 
physics. Alternatively, one could also imagine to perform systematic analysis of the 
differential Galois group of the corresponding linear differential operators in order to 
obtain a deeper understanding of these operators. To some extent differential Galois 
group analysis and p-curvature calculations (see below) are very close. In pratice, the 
analysis of the differential Galois group requires much more time and becomes very 
difficult to perform for linear differential operators of order larger than four, which is 
the vast majority of the operators we encounter in physics. In contrast, p-curvature 

§ Generically the indicial polynomials of Fuchsian ODEs, with coefficients in C[a;] , do not have integer 
coefficients. In our examples originating from a lattice problem we do have integer coefficients (see 
[Appendix C\ . 

jj In this arithmetic framework, this can be seen as a consequence of the Kronecker's theorem. When 
the polynomial coefficients of Fuchsian ODEs have integer coefficients the critical exponents are 
algebraic numbers. In a globally nilpotent context the critical exponents necessarily reduce to integers 
modulo every prime. Algebraic integers reducing to integers modulo every prime are necessarily 
rational numbers (Kronecker's theorem |34l 1351 |36| V 
% Le bourgeois gentilhomme (Moliere). 

f Which is not a surprise for Yang-Baxter specialists, see for instance [39) . 
ft In fact, this operator was first introduced by the Chudnovskys. 
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calculations provide more partial information (only a finite amount of primes can be 
checked), but are easy to perform, simple, and effective. 

We will consider many n-fold integrals, or generating functions, that originate 
from the square Ising model or from enumerative combinatorics. Some quantities are 
not naturally expressed as n-fold integrals of an algebraic integrand: the discovery of 
the global nilpotence of the corresponding (minimal order Fuchsian linear differential) 
operator has to be seen as a strong indication that they are DFG (can be expressed 
as n-fold integrals of an algebraic integrand). Other quantities, like noticeably the 
X^"^ contributions of the magnetic susceptibility of the square Ising model [21 [1], are 
defined as n-fold integrals of an algebraic integran 41. The integrand of the algebraic 
functior lffi can be chosen continuous and single-valued on the torus of integration: 
they are indeed |44j a, family of periods (in other words DFG). 

In this last case, our purpose is not to give another proof of this global nilpotenc^, 
but to see how these linear differential operators manage to be globally nilpotent. The 
corresponding minimal order differential operators can be of a quite high order (for 
example order 33 for x^^\ see [5]), but are always factorisecU into linear differential 
operators of smaller orders. Necessarily, all these factors have to be globally nilpotenl[i|. 
The global nilpotence of order-one linear differential operators is easy to see: their 
wronskians are iV-th roots of rational functions. The global nilpotence of order-two 
linear differential operators is much more interesting: are they "Dworkian" operators 
(see before), and, more specifically, do they correspond to 2F1 functions, or do they 
correspond to counterexamples similar to Krammer's, where Heur0 functions occur, or 
to the more general counterexamples of Dettweiler and Reiter [42]? We will also display 
many globally nilpotent linear differential operators of order three, four, etc. Are they 
reducible to the global nilpotence of the previous order-two operators, because they 
are equivalenlj^, to symmetric squares, symmetric cubes, ... of globally nilpotent linear 
differential operators of order two? Do they correspond to selected „+iF„? We will 
see that the answers to these questions are quite non-trivial, and shed an interesting 
light on the very nature of the globally nilpotent operators emerging from physics. 

From an "experimental mathematics" viewpoint, checking the global nilpotence 
of linear differential operators amounts to studying these operators modulo as many 
primes as possible |48l [49], more precisely by calculating the p-curvature of these 
differential linear operators mod prime for different primes (see below) . Along this line 
it is worth recalling that, in a previous paper [5], we performed massive calculations on 
series expansions of many n-fold integrals. These massive calculations were performed 
modulo various prime numbers and enabled to get many highly non-trivial exact 



§ The corresponding linear differential operators are thus holonomic [21 1431 . because the integral of 
a holonomic D-module is necessarily holonomic. 

ft In the integration-variables e^''^'^^ and in the "parameter" w, see |[T0)| below. 

f The n-fold integrals (over a closed n-cycle) of rational expressions are necessarily DFG. More 
generally, n-fold integrals (over a closed n-cycle) of rational expressions on some algebraic 
variety [II 1451 146] are necessarily DFG. 
f As direct sum factorisations or straight factorisations. 

II The product of globally nilpotent operators is necessarily globally nilpotent. More precisely, 
the characteristic polynomial of the p-curvature of the product operator is the product of the 
characteristic polynomial of the p-curvatures of the factors (see Theorem 5. in 1321 or Corol. 2.1.3. 
in [40]). 

^ They are straight generalisations of hypergeometric ODEs, four singular points [0, 1, a, 00] replacing 
the three points [0, 1, 00] of 2^1 (see the Heun functions generalizing 2^1)- 

tl In the sense of the equivalence of linear differential operators i47j . We refer to this (classical) notion 
of equivalence of linear differential operators everywhere in this paper. 
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results on these physical quantities, thus showing that modulo prime calculations 
are not artificial or academic: they are actually a very powerful, and efficient, tool 
to get highly non-trivial exact results and they are possibly the only way to get 
some "extreme" results in physics. The p-curvature calculations performed here 
are a natural extension of the series and ODE mod prime calculations performed 
in [5]. Let us recall that the linear differential operators that annihilate our n- 
fold integrals, factorise in operators of much smaller order. In this paper, we will 
systematically calculate the p-curvature of these linear differential operators (for 
moderate p) but also of each of the differential factors in their factorisation (in direct 
sums and in products of differential operators). Except in the examples for which 
we do not have an n-fold integral representation of the holonomic function, these 
two set of p-curvature calculations are not performed to check a global nilpotence 
that we know to be a simple consequence of the integral (of an algebraic integrand) 
representation, but to get a global understanding of these ODEs beyond the usual 
local analysis (singularities, exponents, formal series expansions, see the notion of 
"accesso ry p arameters" [50l [5T1 [52l [53] below) and get more precise details on these 
operator af ji . 



The paper is organised as follows: we first recall a few n-fold integrals and 
some basic facts on global nilpotence. The factorization of the linear differential 
operators annihilating these holonomic n-fold integrals will provide a bunch of non- 
trivial examples of globally nilpotent operators of growing orders. We successively 
consider the global nilpotence of such operators of order two, three and four. This 
study will provide a deep understanding of these global nilpotence from the discovery 
of the underlying structures of selected algebraic varieties (hypergeometric functions 
with a Hauptmodul pull-back, various modular structures, etc.). We will finally show 
that there is clearly a generalization of global nilpotence to be discovered in the scaling 
limit of our lattice models. We will conclude with a systematic program of analysis of 
n-fold integrals in theoretical physics. 

In our learn-by-examples approach the variable in the is generally called 

X, except when we need to recall previous results on the Ising model, where the 
variables were called w or t, or when we need to introduce some change of variables. 

2. Recalls on Puchsianity and global nilpotence. 

2.1. Recalls on Fuchsianity for lattice problems. 

When a (minimal order) linear differential operator with polynomial coefficients 
is discovered for series expansions in lattice statistical mechanics, or enumerative 
combinatorics on a lattice, one always finds out that it is a Fuchsian differential 
operator. This lattice property is not true in the scaling limit (see section pop below). 
The regularity of two singular points of the ODE, and cx), is a simple consequence of 
the well-known existence of various kinds of series expansions (low-temperature, high- 
temperature, high-field, large q expansions, . . . ). The fact that the other singularities 
are regular may seem more mysterious at first, except if one remarks that these series 

ft Like the characteristic and minimal polynomial of the p-curvature, the Jordan-block reduction of 
the p-curvature, hopefully in order to get some hint on the factorisation or direct-sum decompositions 
of these operators. 

^ To avoid multiplying the notations, we will sometimes use the same notations for different operators 
when there is no possibillity of confusion. 
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with rational number coefRcients, or even integer coefficients, liave a finite radius of 
convergence, and, in fact, are G-functions (see section ([7]) below). 

It is important to remark that the Fuchsian ODEs one encounters in lattice 
statistical mechanics, or enumerative combinatorics are not the most generic Fuchsian 
ODEs, but very selected ones. One inherits from their lattice origin the fact that their 
polynomial coefficients have integer coefficients, that all the indicial polynomials of all 
the singular points have integer coefficients, and, thus, that all the critical exponents 
are algebraic numbers, but not necessarily rational ones (see [Appendix C] ). 

Let us consider such an order q Fuchsian linear operator. Denote Xk the n regular 
singularities, including the apparent ones and excluding the point at infinity, and p^^^ 
the local exponents corresponding to the singularity Xk- It is straithforward (using 
Fuchs' relations) to write the rational coefficient in front of the {q — l)-derivative in 
terms of the local exponents of the various singularities: 

^ + ^ x-xu + . - dx 

k=l 

The local exponents at the regular singular point Xk are roots of an indicial equation 

(k) (k) (k) _i 

which is a polynomial in p with integer coefficients, Oq + a\ p + ■ ■ ■ + agj_-^ p''^ + 

a^q^p'^ — 0, and, thus, the sum X]j=i Pfe^ = ~^^q-il^^<t\ associated to the singular 
point Xk, is necessarily a rational number. One sees, as a consequence, that there 
exists an integer N such that the N-th power of the wronskian 



W{x) = n ~ = ' pi^'-9(9-l)/2^ (1) 

k 



is a rational function. 

All the examples of Fuchsian ODEs displayed in this paper (see in particular (PT|) 
below) have wronskians that are iV-th roots of rational functions, and, as we just saw 
it, this is straightforwardly inherited from the underlying lattice. Lattice statistical 
mechanics and enumerative combinatorics naturally provide Fuchsian ODEs with N- 
th roots of rational function wronskians, and algebraic numbers critical exponents. 

We will show, in the following, that an even more selected set of Fuchsian linear 
differential operators naturally occurs in theoretical physics, the globally nilpotent 
operators: the previous, at first sight, algebraic numbers critical exponents have 
necessarily to be rational critical cxponentfj 



2.2. Recalls on global nilpotence 

A linear differential homogeneous equation of order q, with polynomial coefficients in 
Q[a;], can always be written as a first order system of homogeneous linear differential 

^ Stricto sensu the rationality of all the critical exponents of a Fuchsian ODE is not sufficient to have 
the global nilpotence property, see 1541 . Global nilpotence is stronger than the rationality of all the 
critical exponents for our "lattice" Fuchsian ODEs. One has conditions on the so-called "accessory 
parameters" 1501 [STI [5^ [5^ . 
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equations : 



Y' = A-Y, Y = 



y 
y' 



or: (^-^)-F=0, (2) 



where the entries of the matrix A are rational functions of x. Instead of studying the 
connection d/dx — A, one can, alternatively, consider for almost anj||| prime p, its 
p-iterate modulo p: 

ipp = (^-^ - A^ , mod p. (3) 

This ijjp is called the p-curvature, and it turns out that this p-curvature for any prime 
number p, is a ¥p{x) -linear operator, so that it can be represented by a matrix whose 
entries are rational functions of x, rather than a linear differential operator. The 
differential system Y' = A-Y yields for the i-ih derivative of Y: 

y(^) = A,-Y, with: A,+i - ^ + A-A, Ai ^ A. (4) 

Katz shows |59| that computing the p-curvature amounts to calculating Ap modulo p 
from the Lie-sequence @ . This can be done by performing p products of qxq matrices 
whose entries are rational functions in ¥p{x) (i.e. rational functions with coefficients 
in Fp where Fp = Z/pZ). These p-curvatures were introduced in the framework of 
the Grothendieck conjecture, to provide "algebraic criteria" for the monodromy group 
of to be finited. 

In the case of order-one linear differential operators, the Grothendieck's conjecture 
was proved by Honda [54]. The fact that the exponents of the various regular 
singularities are rational numbers can be seen as a consequence of the Kronecker's 
theorem which says that any algebraic number which reduces to integers modulo almost 
every prime is necessarly a rational number [341 1351 The conjecture can also be 
proved in some particular cases: by Dwork for ordinary hypergeometric equations |40) . 
and by Katz for Gauss-Manin differential equations (see [47[ [48] for more details on 
the second order linear differential equations when one has only three regular singular 
points, like 0, 1, oo for hypergeometric functions). The conjecture is still open for 
general second order operators. 

Rather than the Grothendieck-Katz p-curvature conjecture, one can consider 
various theorems by Katz [55], in particular Proposition 9.3 in fbW , which shows that 
the reductions modulo p of the Lie algebra of the differential Galois group contain the 
p-curvatures Vp- 

Beyond the situation of the Grothendieck's conjecture where these p-curvatures 
vanish, another highly selected situation corresponds to the case where these p- 
curvatures are nilpotent modulo p, for almost all primes p (for all primes except 
a finite set o/ primes). In that case, the linear differential operator is called globally 

I Almost all, here, and in the following, means for all primes except a finite set o/ primes. 
§ The Grothendieck-Katz p-curvature conjecture is a problem on linear ordinary differential 
equations, related to differential Galois theory. It is a conjecture of A. Grothendieck from the late 
1960s, and apparently not published by him in any form; it has been publicised, reformulated and in 
some cases related to deformation theory proved by N. Katz in a series of papers 1551 1561 l57l 1581 . The 
question is to give an arithmetic criterion for when there is a full set of algebraic function solutions. 
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nilpotent. A globally nilpotent linear differential operator is necessarily a Fuchsian 
linear differential operator, but it has many more strong remarkable structures. For 
instance, all the exponents of all its various regular singularities are rational, but the 
reciprocal statement is not true: a Fuchsian linear differential operator with rational 
exponents is not necessarily globally nilpotent. Global nilpotence is a stronger structure 
than having regular singularities with rational exponents [54j . It is a very strong 
arithmetic property with a large number of remarkable consequences: for instance 
modulo any prime p the Fuchsian linear differential operator factorizes, and for almost 
all primes, it factorizes into linear differential operators of order one, each operator of 
order one having rational solutions modulo p. Such a property is quite well illustrate43 
in Appendix H of [601161] on n-fold integrals related to Apery's analysis of C(3)- In 
that case, we even have a factorisation into order-one linear differential operators on 
the rationals Q and not only modulo (almost all) primes. The fact that the solutions 
of these order-one linear differential operators are actually rational solutions modulo 
primes, is clearly reminiscent of the occurrence of wronskians that are A'^-th roots of 
rational expressions. 

Global nilpotence is often said to suggest a "deep geometrical interpretation", 
namely that the solutions of a globally nilpotent linear differential operator can be 
interpreted as periods of some (hidden . . . ) algebraic variety, suggesting more or less 
a Gauss-Manin connection [7] interpretation for these linear differential operators. 

Beyond the linear differential operators associated with Apery's analysis of C(3), 
almost all examples of globally nilpotent linear differential operators correspond to 
hypergeometric functions, and other Katz's rigid local systems [62], for which an 
interpretation of the solutions as periods of an algebraic variety plays a central role. 
Within the known examples, the overlap between hypergeometric functions (and their 
simple generalisations) and global nilpotence was so large that Dwork proposed a 
conjecture [40] that all the globally nilpotent linear differential operators correspond 
to hypergeometric functions up to simple transformations. This conjecture was ruled 
out by Krammer [41]. Therefore, at the present moment, beyond the fact that it is a 
highly remarkable arithmetico-geometric selected property, one can say that one does 
not have a complete understanding of global nilpotence. 

In the following we are going to find globally nilpotent linear differential operators 
corresponding to various n-fold integrals that occur naturally in the case of the off- 
critical lattice Ising model, or corresponding to enumerative combinatorics for which 
no n-fold integral representation is yet known. We will also explore situations that are 
precious to understand namely other n-fold integrals that naturally occur in particle 
physics (Feynman diagrams [63]), like, for instance, some selected scaling limits. 

2.3. Krammer' s counterexample. 

Let us recall briefly Krammer's counter-exampl Jfjl to Dwork's conjecture [41] |42] 
which comes from the periods of a family of abelian surfaces over a Shimura curve 



pi 



\ {0,1,81, oo}: 












1/2 -1/2 



^ Note a misprint in 1601 one should read In Ai , instead of Ai , in the equations defining the Ai after 
equation (H.2) in leOlleTI . 

It The uniformizing hnear diff'erential equation of an arithmetic Fuchsian lattice |41) . 
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4/9 -1/2 



^81 



1 

1/2 



yielding the second order operator (on the first component of the vector Y): 



Oi 



1 1 
- H 

X X — 



1 



81 



(5) 



(6) 



9 



ax 



(7) 



18 {x - 81) {x - 1) • a;' 

or the second order operator (on the second component of the vector Y) 

O2 = 18a; • (x - 1) (a; - 9) {x - 81)^ • DI 

+ 27{x- 81) • (a;^ - 123 x^ + 1491 x - 729) • 

+ {x^ + 549 x^ + 13203 a; - 1003833) . 

The two hnear differential operators Oi and O2 are, of course, equivalent and the 
squares of their wronskians are simple rational functions. We have calculated their 
p-curvatures and confirmed that they are globally nilpotent. 

The general solution of Oi reads in terms of Heun functions [M] : 

H ■ Heun{81, 1/2; 1/6, 1/3, 1/2, 1/2; x) 

+ A • x^/^ . Heun{81, 21; 2/3,5/6, 3/2, 1/2; a;), 

where fu. and A are two constants. The differential Galois group of ^ (or ([U, ([7])) 
is a central extension of SL{2, C). Calculating the indicial polynomials of Oi at the 
various singularities, one finds the indicial polynomial (6r— l)-(3r— 1) for t — 00, 
and r • (2r — 1) for the singularities t — 0, 1, 81. These Heun functions cannot be 
reduced |41] to hypergeometric functions 2-F1 (up to multiplication and some pull- 
back). Along this line, it is important to note that, generically, a Heun function does 
not correspond to a globally nilpotent second order differential operator. For instance 
we calculated the p-curvature of a second order operator very similar to ([5]): 



1 



1 



81 



81 - 28 a; 



2 (x- 81) (a; - l)x' 



a; — 1 X - 
which has as solutions the Heun functions 

/X • Heun{8l, -81/2, -7/2, 4, 1/2, 1/2; a;) 

+ A • • Heun{8l, -20, 9/2, -3, 3/2, 1/2; x) 
and we found that ([5]) is not globally nilpotent. 



(8) 



3. Global nilpotence of a few n-fold integrals of the Ising class. 
3.1. Recalls of a few n-fold integrals of the Ising class. 

The susceptibility of the Ising model can be written [SS] as an infinite sum of n- 
fold integrals. These n-particle contributions x^"-* are given by (n — l)-dimensional 
integrals [66l [671 EHl that read 

, n— 1 „27r 11 " 
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wher^ 



and 



with 



n h h - 2sin((0i - <j)j)/2) ■ ^Jx^x] 

_L tXJ tJU n 



(10) 



1 < i < J < n 



= ]^?,h (11) 



i-m 



2u; 



) + y'(l-2wcos(0,))^-4u;2 



(12) 



2wcos( 

2u; 



^0, = O (13) 
(1 - 2wcos(0i))^ - 4^2 j=i 

valid for small w and, elsewhere, by analytical continuation. We actually found [31 [H 
[551 [70] the linear ODEs for some of these holonomic n-particle contributions namely 
X^^^ I X^^^ and, modulo a prime, for x'^' . From a arithmetic Gevrey series and G- 
function viewpoint it is worth noticing that the series expansion of the x'""\ the 
variable w, are series expansions with integer coefficients: 

X^^^w) = 2" • w"' • (l + An^-w^ + 2- (4^4 +13n^ + 1) ■ 

+ ^•(n2+4)(4n'* +23n2 + 3)-w^ + •••), (14) 

where the w'^ coefficient is valid for n > 3, the coefficient is valid for n > 5, 
and the coefficient is valid for n > 7. Note that the coefficient is always an 
intege^ 

In previous publications [5D1 |7T], we also introduced some integrals of the so- 
called "Ising class"[3. We considered several kinds of integral representations (one- 
dimensional and multidimensional) of these holonomic functions which belong to the 
Ising class [72] . Again we obtained the linear ODEs of these sets of integrals for the 
first values of n, through series expansions ^Tl|60]. In [60] a detailed analysis of the 

(n) 

multiple integrals was performed. These rt-fold integrals amount to removing 
the Fermionic factor G^"^ in ([9]), so that one introduces: 

Even simpler integrals (over a single variable), were also introduced and 

(n) 
D 



denoted [H] 



(„) , , 1 2 f^^ dtp 1 ,^„^ 

D ^""l = 2^ 1 -a;"-i(0)-a;((n-l)</>)' ^ ' 

where x{(l)) is given by 

f The fermionic term g'"' has several representations |67| . 

X It would be interesting to get much longer series expansion like I I14I I. valid for arbitrary n, to see 
if these successive rational functions of n are actually functions of ri^ . 

tl The terminology integral of the Ising class has been proposed by Bailey, Borwein and Crandall 
in [TD. 
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3.2. Results on global nilpotence of these n-fold integrals of the Ising class. 



Let us display here our results for the calculations of p-curvatures for the minimal 
order ODEs of n-fold integrals p^ . pB]) of the Ising class. 

We have calculated (modulo the first thousand primes) the p-curvature of the 
order-six linear differential operator Lg occurring in x^^"* (see ([22l) in section (|4.1 



and [70]), as well as the Jordan block reduction of the 6x6 p-curvature matrix, and 
found that the characteristic polynomial of the p-curvature reads T^. This 6x6 
Jordan block reduction can be compared with the two 3x3 Jordan block reductions 
corresponding to (the p-curvature of ) an order-three operator Z2-Ni that right-divides 
Lq, and another order-three operator that left-divides Lg (see section (|4.ip below). 
They read respectively (in block form): 



A 





B 



A 



' 





" 




" 


1 


' 








1 


, B - 








1 























The explicit Jordan-block form of the p-curvature is quite reminiscent of the 
factorisation of the operator Lg (see section (|4.1I) below) in an order-three linear 
differential operator, and another order-three operator, itself product of an order-two 
operator and an order-one operator. Note, however, that one should not extrapolate 
beyond simple product factorisations: the Jordan-block form of the p-curvature gives 
systematically a misleading prejudice of direct-sum structures that do not exist. 

The global nilpotence of the order-ten Fuchsian linear differential operator [69] 
for x^^\ is confirmed by the calculation of the p-curvature for all the primes up to 
p < 809. The p-curvature has been found to be nilpotent for all these primes. 

The global nilpotence of the order 5 and 6 Fuchsian linear differential operators 
for (f^' and <i>^'' is confirmed by the calculation of the p-curvature: we have calculated 
the p-curvature for all the primes up to p < 809 and it has been found to be 
nilpotent for all these primes. The characteristic polynomial of the p-curvature of 
the (globally nilpotent) linear differential operator of has been found to be 
(its minimal polynomial being T^). The characteristic polynomial of the p-curvature 
of the (globally nilpotent) linear differential operator of has been found to be 
(its minimal polynomial being T"*). For the calculations are drastically larger, 
but, from a probabilistic algorithm, we found that the characteristic polynomial of the 
p-curvature of the linear differential operator of is T^^. 

The minimal polynomial of p-curvature of ^^^^ and is T*. The characteristic 



polynomial of the p-curvature of and $^ is T^. The minimal polynomial of 

p-curvature of <i>2^ is T^. Recall that the characteristic polynomial of the p-curvature 
of a globally nilpotent operator of minimal order N equals . 

3.3. Other n-fold integrals of the Ising class. 

Other n-fold integrals (corresponding to the susceptibility of a square Ising model for 
which a magnetic field is located only on spins on a particular diagonal of the square 
lattice) were introduced in ^73) . For instance, for T < T^, they read: 

Jo Jo fc^i 



(n!) 



1 + t" Xi • • • X2n 
1 - t" Xi • • • X2n 
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X 



n 



a;2j-l(l - X2j)il - tX2j) 



-I, 



1/2 



7=1 Va;2i(l - a;2j-i)(l -tx2j- 
^ n n (1 -*a:2j-ia;2fc)"' (17) 

l<j<n l<k<n 
l<j<k<n 

and another similar formula for x'"^^^^\t). They are holonomic functions and their 
corresponding Fuchsian linear differential operators were given in [73] . Again the 
calculations of the p-curvature of the corresponding linear differential equations of 
minimal order for Xdl confirmed their global nilpotence (see (j6.ip below 

and [Appendix D| . 



3.4- DDEs for two-point correlation functions and form factors. 

Many simple linear ODEs of various orders were obtained for the two-poinl[^ correlation 
functions of the square Ising model [74j . The two-point correlation functions were 
found to be polynomials (with rational function coefficients) of complete elliptic 
integrals of the first and second kinds: their global nilpotence is, thus, a straight 
consequence of their hypergeometric nature. 

Along this correlation function line, we can recall the linear differential operators 
Fj{N) we obtained for the form fac^or Jffl /j^''jv '^^ (off-critical) square Ising 
model |75j and, in particular, their russian-doU structure. 

The linear differential operators i^2n+i(A^), which annihilate the form factors 
^n^n'^^ have a "russian-doU" structure. They are such that: 
^^i(iV) - L2{N). 

F^iN) = L4{N)-L2{N), (18) 
F5{N) = Le{N) ■ L4N) ■ L2{N), ••• 
where the differential operators Lr{N) are of order r. The first one reads: 

L2{N) = Dj + ' .Dt - —+ , . - (19) 

^ ' * (t- l)t At Ait-l) At^' ^ ' 

and the expressions of L4{N), Lq{N), Ls{N) and Lio{N) are given in [75] . 

Thus we see that the linear differential operator for /^"^ rightdivides the 

differential operator for fl^"j^^\ n < 3. Similar relations occur for the F2„(-/V)'s. 
We conjectured [75] that this property holds for all values of n. We thus have a 
"russian-doU" (telescopic) structure of these successive linear differential operators. 

Again, these form factors were found to be polynomial (with rational function 
coefficients) of complete elliptic integrals of the first and second kinds: the global 
nilpotence of the corresponding operators is, again, a straight consequence of their 
hypergeometric nature. 

§ Or could have been obtained for any Af-point correlation functions. 

ft Coefficients in of C{N, N; A), the A-extension |75| of the two-point correlation function 
C{N, N). 
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3. 5. Modular ODEs for lattice form factors. 

Along this correlation function line, it is also worth recalling the Fuchsian linear ODEs 
we found [75] for some A-extensions C{N, N; A) of two-point correlation functions 
of the (off-critical) lattice Ising model for selected values of the parameter A. As 
examples of these Fuchsian linear differential operators, we found, for instance, that 
C-(-/V, A^; cos(7r/4)), for A^ = 0, 1, 2, • • •, are annihilated, respectively, b-\ lffl 

^^1= it-lft-D^ +lit-l)i3t-2)-Dt 

it-irt.D^ + ^i^i}^^.Dt-ii+^, 

1 V y t T g 256 16' 

7 



4^/"! = (t- 8){t -ift-Dj +--{t- Vfii"- -2t+ 16) • Dt 

8 

209 25t 1 



256 16 2' 



(20) 



One has homomorphisms between the linear differential operators with same parity: 
the Ln^^^ for n odd (resp. even) are equivalent [1H|- We also have higher order 
ODEs like, for instance, the order- four linear differential operators [75], denoted ij^^^' , 
corresponding to C_ (A^, A^; cos(7r/3)). 

We have calculated the p-curvatures of all these irreducible linear differential 
operators and seen that they have zero p- curvatures. Not surprinsingly the 
corresponding wronskians associated to these A-extensions of two-point correlation 
functions are A^-th root of rational functions and read: 

W^(4/'l) = (l-t)-3/8.t-3/4, W^(if/"1) = (l-t)-3t-2, 

W^(4l/4]) ^ (t-8)-8/7.(l_t)-15/7.i-2/7^ ... (21) 

W^(4/3]) ^ _ . ^-11/3^ T4^(i[l/31) = t-14/3 . _ ^)-8/3_ 

t^(4/3]) ^ (11 + 21 i) • (1 - tf/^ ■ 

The fact that the 8-th power (instead of the square in most of the examples of 
this paper) of W{Lq^^^), or the 7-th power of 14^(^2^^'''), is rational is in agreement 
with the interpretation of A we gave in [7S]. These Fuchsian linear ODEs actually 
correspond to algebraic function^, and are often, explicitely, associated with modular 
curves. 



3.6. More zero p-curvatures: Joyce's Fuchsian ODEs. 

Finally, along this zero p-curvature line, it is also worth recalling the large set of 
higher order Fuchsian ODEs obtained by Joyce [7i jTZljlHl [71 [Ml [51] . There are not 
so many examples of Fuchsian linear ODEs of high order in the litterature. Joyce has 
been one of the few authors to provide such non-trivial examples. We have calculated 
the p-curvature of a large set of these Fuchsian linear differential operators, namely 
(42) in [77], (85-86) in [TS], (5.22) of [75], (2.16) of [73], . . . We found that they are 
more than globally nilpotent: their p-curvature is zero for almost every prime (for all 
primes except a finite set of primes). They have a basis of algebraic solutions that can 

ft Note two misprints in |75| for Lj^^*^ and Lj^^^l^ corresponding to the Dt coefficient. 

f Note that maple also solves these ODEs in terms of (algebraic) hypergeometric or Heun functions 

(see 1I35J, (I37J below). 
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be expressed in terms of simple Legendre-P functions and simple algebraic functions. 
These Fuchsian linear ODEs actually correspond to algebraic functions, and are often, 
explicitely, associated with modular curves. 



4. Global nilpotence from the global nilpotence of the factors. 
4-1. Revisiting the global nilpotence of x^^^ 

The minimal linear differential operator for x'^'^^ is an order-seven operator Ly which 
can be written as the direct sum of the order-one linear differential operator for 
X^^"^ and an order-six linear differential operator for 2x^^'^ — X^^\ namely Lq which 
factorizes into an order-three linear differential operator Y3, an order- two linear 
differential operator Z2, and an order-one linear differential operator A^i: 

Le = Y^-Z2-Ni. (22) 

The explicit expressions of and Z2 are given in Appendix A of [7D] . 

We have the following, almost obvious remarli3|, that the product of globally 
nilpotent operators is necessarily globally nilpotent. The global nilpotence of the order- 
one linear differential operator A'^i is obvious. Furthermore, we found the remarkable 
result that the solutions of Y3 are quadratic expressions of the complete elliptic 
integrals of the first and second kind K and E (see appendix B of [7D]). From 
a differential algebra viewpoint this amounts to saying that Ya is equivalent to the 
symmetric square of the second order linear differential operator Le corresponding 
to E (see [HI])- Since hypergeometric functions correspond to globally nilpotent 
operators, Ya is therefore globally nilpotent. 

The global nilpotence of the linear differential operator for x*^"^^ thus reduces to 
the global nilpotence of the second order linear differential operator Z2. The linear 
differential operator Z2, is an example of globally nilpotent operator which does not 
straightforwardly reduce to hypergeometric functions up to change of variables (pull- 
back) and multiplications. 

In fact, a simple right-multiplication of Z2 by h{w) = 1/(1-1- Aw)/{1 — 4u')^, 
enables one to get rid of the singularity w = —1/4. Instead of the solutions F{w) of 
Z2, this amounts to considering the second order linear differential operator Z2 with 
solutions F{w)/h{w). Denoting = d/dw, this linear differential operator reads: 

Z2-h{w)^ q2 ■ Qapp ■ D1 + qi ■ D^, + 24:W ■ go, 



with: 



h{w) 



(72= w(l-w)(l-4w) (l + 2w)(l + 3w + 4w2), (23) 



90 



l-2w + w^ ~216w^ ~ 336 w'^ + 1656 + 1040 
- 2560 - 6400 - 6144 w^, 

1 - 7 w - iw^ - 47 + 36w^ + 280w'^ + 160 + 256 u;^, 
qapp ^ 1 -3w -18w^ +10Aw^ +96w^. 

The polynomial qapp corresponds to apparent singularities. All the other singularities 
are regular singularities remarkably with integer exponents and all yielding logarithmic 

^ If one takes as a definition of global nilpotence the factorization modulo (almost all) primes of the 
operator in order-one operators with rational function solutions modulo primes (see |32] and Lemma 
0.6.2 in [40]'). 
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behaviours. This can simply be seen from their corresponding indicial polynomials 
and the formal series around these singularities. The differential Galois group of Z2 
is SL{2, C), a consequence of the existence of logs in the formal solutions around 
singular points, together with a wronskian being rational, and the operator being 
irreducible [83] , 

Calculating the indicial polynomials at the various singularities, one finds the 
following critical exponents: 

p=0, 0, for: w^O, 1/4, ~^ t 

8 

p= 0,2, for: w = 1, -1/2, 

p = -1, -2 for: w = 00. (24) 

If one wants to see the solution as hypergeometric functions up to a pull-back, the 
change of variables to be done must try to "wrap" all these seven singularities onto 
only three ones: 0, 1, 00. Naively one can thing of wrapping the singularities 
according to the previous roots of indicial polynomials (|24p . namely w — 0, w = 
1/4, (-3 ±iV7)/8 ^ 0, w = 1, w = -1/2,^ 1/4, w = 00 ^ 00, however, 
since the critical exponents are all integers, these seven singularities have, in fact, to 
be considered on the same footing. This will be discussed in more details in section 
(j5.ip . where the explicit solution of Z2 will be given in term of a modular form of 
weight one. 

4.2. Revisiting the global nilpotence of x*-^-* 

For the order-ten linear differential operator of x*"^"* have similar calculations. The 
order-ten linear differential operator Liq of x^*-* is the direct sum of the order-two 
linear differential operator for x^^^ ^'Hd of an order-eight linear differential operator 
which factorizes [69l EO] into an order-four operator and four order-one operators (see 

(F.4) in feg;): 

Cs = M2 ■ L25 ■ L12 ■ ■ Lf), (25) 

The global nilpotence of the order-one operators L25, L12, L3 and Lq is a simple 
consequence of the fact that all these operators are of the form + R' /R, where 
R denotes a rational function, or the square of a rational function (simply related to 
the wronskian of these operators), and R' its first derivative with respect to x. These 
R functions read respectively for L25, Lq, L12 and L3: 

R25 = ^ ^— ^ ^, Rq = 1, (26) 

{U-12x + 9x'^ -5x^)-x 2 + 3x'^+x^ 

R12 — —7Z :^ — T. TTTTZ s^/^ , R3 



(2 + 3a;2+a;3)(l-x)V2 ' ' (l-x^-x^/'^' 

The order-two linear differential operators Ni, i — 1, • • • , 9, sketched in [70] and 
which happen in other factorizations of Liq, are all equivalent to the order- two linear 
differential operator [70] Nq 

Nq^ Dl - + -r- 1 ' 27 

1— X X 4a; 1 — a; 

having hypergeometric solutions (corresponding to x'^'')- 

X(2) = . 2Fi([3/2,5/2],[3],x). (28) 
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Most of the "complexity" of x*'"'"' is thus "encapsulated" in the order-four linear 
differential operators M2 (or equivalently Mi of [70]V 

In a remark (in Appendix B of [7^ page 27), we mentioned the fact that the 
solutions of M2 can be expressed as linear combinations of products of complete 
elliptic integrals. One can make this statement more precise. Let us introduce the 
symmetric cube of the linear operator |70j Nq (associated to x*"^"*)! order-four 
linear differential operators M2 (or equivalently Mi) is equivalent to the symmetric 
cube Sym'^{No). There exists two order-three interwinners Ii and I2 such that: 

h-Sym^No) = M2 ■ h- (29) 

In other words, the solutions of M2 are cubic (homogeneous) polynomials [7D] 
of the two solutions of Nq (hypergeometric functions). Therefore M2 is globally 
nilpotent, and consequently the order-ten linear differential operator for x^^^ is also 
globally nilpotent. One sees that, paradoxically, the global nilpotence of the linear 
differential operator for x''^^ is much simpler to understand than the global nilpotence 
of the linear differential operator for x''^': which is a consequence of the, after first 
sight, subtle global nilpotence of Z2. 

4-3. Revisiting the global nilpotence of $^ . 

In order to revisit the global nilpotence of (see ([T5|)'). let us study the factorisation 
of the corresponding order- five Fuchsian linear operator given in [60| . This order- five 
Fuchsian linear operator factorizes into an order-three and an order-two Fuchsian 
linear operators: 

L^(3, = M3 • M2, (30) 

where the second order operator AI2 is given in [Appendix B.2[ The wronskian of the 
order-two operator AI2 is a simple rational function. The indicial polynomials for M2 
yield integer critical exponents, when the formal series solutions show logarithms for all 
the singularities, the second order operator M2 being irreducible. Consequently [83], 
the differential Galois group of M2 is SL{2, C). 

We have calculated the p-curvature of M3 and M2 and found that their 
characteristic polynomial equals their minimal polynomials, being respectively 
and for almost all primes. These two linear differential operators are thus globally 
nilpotent. 

More remarkably, we actually found that M2 is equivalent to Z2 of X^'^\ (see 
(j41l) below). We also found that M3 is equivalent to of x^^\ or equivalently 
Sym^ {Q e) , the symmetric square of the linear differential operator Qe corresponding 
to the complete elliptic integral E{4:x): 

Qe= Dl + ^ + ——^^ — . (31) 

X {I - Ax) (1 +4x) 

Therefore x^^'' and have extremely close structures. 

As usual the order- two intertwinners Hi and H2 (or H[ and H2) in the 
equivalence of M3 and Y3 

M3 Hi = H2- Sym^QE), H[ ■ A/3 = Sym'iQE) ■ H'2, 



t Up to the change x = iw. 
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are themselves equivalent and, again, there exist two order-one intertwinners Ki and 
if 2 (resp. K[ and K'^) such that Ki ■ Hi = H2 ■ K2, and K[-H[ = H'^-K'^, which 
are again equivalent. One thus has a "tower" of equivalent differential operators. It is 
important to note that the "intertwinning" operators in such a "tower" of equivalences 
are not necessarily globally nilpotent! 

4.4- Revisiting the global nilpotence of 

The global nilpotence of can be understood from the factorisation of the 

corresponding order-six Fuchsian linear differential operator given in [60] . It factorises 
as follows: 

L^(4) = A'U-Pi-Qu (32) 



where the order four operator A/4 is given in [Appendix B.3 and the two order-one 
operators Pi and Qi read: 

Pi= + i Aln((a;-4)(x-l)2.x2 

Qi = + 1 Ain((x-1)2. 

The global nilpotence of Pi and Qi corresponds to the fact that the square of the 
wronskians of these first order operators are simple rational functions. 

The order-four operator M4 is irreducible but it is actually equivalent to 
Sym^{LE), the symmetric cube of the linear differential operator Le having E{x^^^) 
(complete elliptic integral) as solution: 

Le = Dl +^ +I7r^- (33) 
X 4, [1 — X) X 

The global nilpotence of ^^'^ is clearly a straight consequence of this last 
equivalence with a symmetric cube of a hypergeometric operator. Comparing 
the factors (and their equivalence) for one finds out, after the remarkable 

identification of structure between and x'^'^K that one has exactly the same 

identification of structures for Cj^-* and for x^''''- Again, slightly surprisingly, the 
global nilpotence of <I>^'' (or x*^**') is straightforward to understand compared to the 

global nilpotence of (or x^'^'') which amounts to understanding the more subtle 
global nilpotence of the second order operator Z2. 

(n) 

4.5. Revisiting the global nilpotence of the 's. 

Similar (detailed factorisations) calculations can be performed for ^'"^\ ^'d' ^ 



and defined by (|16p . They are displayed in Appendix E Again one finds out 



that the global nilpotence of many of the factors occurring in the factorisations can be 
explained by the occurrence of complete elliptic integrals of the first or second kind. 
The example of ^^^^ in Appendix E| corresponds to a much more interesting situation. 



also occurring with , where we contemplate in the factorisation of the operator for 

4>^'' , a remarkable combination of a globally nilpotent order-two operator associated 
with complete elliptic integrals of the first or second kind, together with a second order 
operator of zero p- curvature associated with highly non-trivial (genus five) algebraic 
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curve (with a group of quaternion as its differential Galois group). Let us illustrate a 

(5) 

similar structure focusing on the global nilpotence of . 
4-6. Revisiting the global nilpotence of • 

(5) 

The linear differential operator for <i>]j is a Fuchsian linear differential operator of 
order five which factorises in three different ways {N4, N3, N2, Ni are second order 
linear differential operators, Li and L2 first order): 

L^,5) = Li-N2-Ni^ N3-L2-Ni^ N3-Ni-D,. (34) 

This means that it factorises into the direct sum of two order-two linear differential 
operators and the operator Dx'. 

L.(5) = Ad4®Dx, with: M4 = N2 ■ Ni, 

where the second order operator Ni is given in [Appendix B.4| The square of the 
wronskian of iVi is a simple rational function. 

The p-curvature of A^i is equal to zero for any prime > 2, therefore, modulo 
Grothendieck's conjecture[§|, it admits a basis of algebraic solutions. These algebraic 
solutions can be written as hypergeometric functions up to a pull-back by a rational 
function and multiplication by a A^-th root of a rational function. They are of the 
form [SlIHg: 

Ri{x)^/^ ■ H[R2{x)], where: P{H[R2{x)], x) ^ 0, (35) 

where Ri{x) and R2{x) are rational functions, H is an algebraic hypergeometric 
function and P{y, x) = is an algebraic curve of degree 4 and of genus 6: 

{I- X - ix^ + A x^f {l + xf {l + 8x + 2Qx^ + lbx^ + A x^) 
X (1 + 2 a; - 4 a;2) (1 - 3 x + (1 + 2 x) (1 - x) • y'' 
-2{l + x){l~x-ix^ + 'ix^) ■ gi-y^ (36) 
+ (1 +3a; +3a;2 - Ax^f [l + x - Qx^ - x^ + x^f ^ 0, 

where: 

= -3-24a;-20a;2 + 255 2:^-|-484a;^-800a;^- 1729a;^ 

+ 1296 + 2236 a;*^ - 1035 - 1004 a;i° -I- 232 a;" -I- IQOx^'^. 
Note that the vanishing of the y'^ coefficient in (1551) corresponds to singularities of 

This genus six algebraic curve should not be confuseqg with genus three curves 
that will be mentioned at the end of section ( [Appendix F.l[ ). 

Using the maple command kovacicsols one can write one of the solutions in 
terms of hypergeometric functions as follows: 

(i?^g,(a;))'^''. 2i^i(hl/4,l/4],[l/2],X), (37) 

§ Namely that zero p-curvatures yield 1141 1551 1561 1571 a basis of algebraic solutions for the linear 
differential operator. 

tt They can, however, correspond to highly non-trivial relations expressing such genus six algebraic 
solutions as linear combinations of complete elliptic integrals of the third kind with a "characteristic" 
(first argument of the complete elliptic integral) associated with a genus three curve. 
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where the argument in the hypergeometric function reads: 
TV 

X = — , where: (38) 

TV = (x-l) (l + 2a;)(l-3.T + .T2)(l + 2x-4x2) 

X (1 + 8 x + 20 a;2 + 15 + 4 a;-*) (1 + 3 a; + 3 - 4 x^f 
X {l + x-6x^ - x^ +x^)^, 
V= 8 (1 + a;) (1 - 3.x2 + 4a;^) -d^o, 
dio = 1 + 8 a; + 8 - 73 a;^ - 148 + 144 x^ + 325 x^ 

- 152 a;^ - 168 a;^ + 39x^ + 28x1°. (39) 
and the rational function i?^(5) (a;) reads: 

Tl 

^$(5) (a;) = 2-^:7' where: (40) 

n= 16 • (a; + 1)2 (1 - a; - 3a;2 +4a;3)2 . d^Q. 

The hnear differential operator (or equivalently TV3) has solutions in terms 
of the complete elliptic integrals E and K . The p-curvatures of the linear differential 
operator N2 (resp. Ny,) correspond to a globally nilpotent operator. 

(n) 

Towards a geometric interpretation of the • 

The DFG (global nilpotence) structure corresponds to the fact that a holonomic 
function has an interpretation as a period of some algebraic variety. Along this line, it 

in) 

is worth noting that in the case of these tj)^) , some closed exact expressions for these 
integrals C'^' can be obtained which give explicit examples of such an interpretation 
as a period. Actually the simple integrals can all be expressed as sums of complete 
elliptic integrals of the third kind, where the characteristic y = y{w) corresponds to 
some selected, and highly non trivial, algebraic curves (genus zero, three, ten, ...). 
The results are displayed in [Appendix F.l| 

5. The second order operator Z2 and weight-1 modular forms. 

The previous section was dedicated to understand the global nilpotence of linear 
differential operators of quite high orders from the global nilpotence of their factors of 
smaller orders. To sum up all the previous results, the global nilpotence of the factors 
fall in three categories: 

- firstly, a straightforward global nilpotence of order-one factors (which amounts 
to saying that their wronskian is a A^-th root of a rational function), 

- secondly, global nilpotence straightforwardly associated with complete elliptic 
integrals of the first (or second) kind, or equivalently 2^1 hypergeometric functions, 
or global nilpotence corresponding to symmetric powers of the previous second order 
hypergeometric operators, or to zero curvature operators with their solutions being 
a basis of algebraic functions (often modular curves, such algebraic functions can be 
written as pull-backs of 2-F1 hypergeometric functions), 

- finally, a set of operators of order two, three, . . . that do not have a basis of 
algebraic solutions and that we have not been able to immediately reduce to 2-^1 
hypergeometric functions, or products of hypergeometric functions. 
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The second order linear differential operator Z2, occurring in the factorisation 
of the linear operator for x'^-* (or equivalently M2 for ^'^''), is a perfect illustration 
of this last situation. The order two and three operators occurring in the operator 
factorization for the three-choice polygons perimeter generating function (see below) 
are other examples that do not immediately fit in a 2F1 "Dworkian" framework. 

5.1. Reducing Z 2 to Heun functions 

We found that the Z2 of x*^'^\ ^^id its equivalent for namely M2 in (|30p and 

(|B.3p . are homomorphic0 (Ui and U2 are two order-one intertwinners): 

Za-f/i = U2-M2, C/i - + (41) 

(l-2a;)(l + 4a;)(l-4a;) ' ^ ' 

ui = X ■ {1 - Ax) {1 + 3 X + Ax^) ■ + {1 + X + 12 a;^ + 48 .x^), 

On the other hand, let us introduce the order-two Heunian operator whose 
solution is Heun{8/9, 2/3, 1, 1, 1, 1; t): 

^< +(l+7^+9Fb)-"' f"^' 

3t-2 



{9t-8) (i- l)t 

a simple change of variable: 

(1-4.)(1-.) ' (43) 
transforms (|42p into the order-two linear differential operator: 

= Dl (44) 
1 - 10 a; -M9 a;2 - 92 a;3 + 12 a;"* 



(1 -I- 3 a; + 4 a;2) (1 - 2 a;) (1 + 2 a;) (1 - 4 a;) (1 - a;) • X 
(l + 7a; + 4x2)(l-2a;)2 



6 



(l + 3a; + 4a;2) (1 - 4a;)2 (1 - a;) 



2 



X 



This operator ((44|) is just the conjugate of M2 by the multiplication by a simple 
polynomial function h{x) = (1 — a;) (1 — Ax)"^: 

h{x) ■ M2 = ■ h{x). (45) 

Together with (|4T|) . this means that the second order linear differential operator 
Z2, corresponding to x''^\ reduces to a Heun operator given by (l42|) . to be compared 
with Krammer's [41] counterexample ^ to Dwork's conjecture. 

One then gets the (selected) solution for Z2 in terms of Heun functions {w replaced 
here by a;, and t is given by ((43l) ): 

r{x) ■ {(I - 9 a;) (1 - 4 a;) (1 - x) • Hg{t) 

+ 8-{l + ix + Ax^)-x- Hg'{t)^, (46) 

>, ( \ (1 +2a;)2 
where: r[x) = -; -15— — — — and: 

(1 - xY (1 +4a;) (1 - 4a;)2' 

Hg\t) = Hg(t) = i/emi(8/9, 2/3, 1, 1, 1, 1, t). 

^ Up to the change x = 4 to in M2. 
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5.2. Reducing Z2 to weight- 1 modular forms 

Recall |86j that the fundamental weight- 1 modular forn^ hj^ for the modular 
group ro(-/V) for = 6, can be expressed as a simple Heun function, 

Heun{9/8, 3/4, 1, 1, 1, 1, — 1/8), or as a hypergeometric function: 

2^3 



{{t + 6)3 (t3 + 18 i2 + 84 1 + 24)3)1/12 

x.^.f[4^Ml];1728 (^ + 9)2(. + 8)3. 



(47) 



42' 12^'^ " (t + 6)3(i3 + i8i2 + 84i + 24)3 

and that ([T71) is solution of the order-two linear differential operator (obtained from 
(ga) by i ^ -t/9) : 

+(^+i+^)-^- + (^ + 8)(^ + 9).r 

Therefore, after some changes of variables, one can see the (selected) solution of Z2 
as a hypergeometric function (up to a pull-back) corresponding to weight- 1 modular 
forn^ {he in [86 ). To sum-up Tix, given by (|44p . has the following solution: 

S ^ (n-M:,y^'\ 2Fi{[^,^];[l]; M,), where: (49) 
^ _ 1 (l-4x)6(l-2:)6 



1728 (1 -|-3a; -t- 4x2)2 (l-f 2a;)6 ' 

a; - (1 +3 a; +4a;2)2(i -|- 2 x)^ (1 - 4 a;)^ (1 - x)^ 



Ma: = 1728 



(1 -f 7x-|-4x2)3 .p3 

P ^ 1 + 237x + 1455x2-f 4183x3 -f 5820x^^4- 3792x5 64 x^ 

Recalling (|^T|) and (|l5)) . the solution of the operator Z2 in terms of hypergeometric 
functions then corresponds to the action of the intertwinner Ui on the solution S/h{x) 
of M2 . The global nilpotence of Z2 can now be understood from this hypergeometric 
function (up to a modular invariant pull-back) structure. 

5.3. Atkin's modular functions 

The modulus Aix in the argument of 2F1 actually corresponds [SB] to the (genus 
zero) modular curve which amounts to multiplying, or dividing, the ratio of the two 
periods of the elliptic curve by 6: 

Mj, ]') = j) = 0, (50) 

obtained from the elimination of z between: 

(z + 6)^ (z3 -1-18 z2 + 84 z + 24)3 

^ " ^^^^^ = ( ^0^2. ^^.3 ' (51) 

z(z + 9) (z + 8) 

' 23 • 32 \ (z + 12)3 (^3 ^ 252 z2 + 3888 z + 15552)3 

JeC - 



z6 (z-F8)2 (z + 9)3 



t The modular form /ig is also combinatorially significant: the perimeter generating function of the 
three-dimensional staircase polygons 1871 can be expressed in terms of (see section US. Il l below). 
The modular form also occurs in Apery's study of ^(3) (see [Appendix A\ . 

§ The simplest weight-1 modular form is 2-Fl([l/12, 5/12], [1], jf^T^T^ ^(t-)2 j-1/12 „here J 
is the Hauptmodul, rj is the Dedekind eta function, and t is the ratio of periods (see (4.6) 
in [88]). It can also be expressed as a linear combination of 2^i([l/12, 5/12], [1/2], 1 — J) and 
2Fi([7/12, 11/12], [3/2], 1 - J) . (1 - j)i/2. 
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Actually, similarly to if one introduces: 

72 X 12^ 

' ^ (i-.x)(i-4.) ' ^^^^ = m ^''^ 

one finds immediately that in (j49p is nothing but the Hauptmodul M^- 

The singularities of the linear ODE of x^^^ (or $^'') correspond to the 
singularities of j{z): 

z{z + 9){z + 8) = (53) 
x-(l + 3x + 4x^)(l+ 2x)^ 
(l-x)3(l-4x)3 
Similarly the "Atkin-dual" [SS] change of variables (see [Appendix G ) : 



,^72^ il-x)il^Ax)^ 

Z X 

gives: 

z'(z' + 9)(z' + 8) = (55) 
(l-a;)(l-4a:)(l + 2a;)^(l + 3x +4a;^) 
x^ 

Recalling the modular Atkin-polynomia^ [52] for ||5D|) (see [Appendix G[ ): 

z- (z + 9)2(z + 8)3 • j - (2 + 6)3(^3 + 18 2^ + 84 z + 24)^ (56) 

one finds out that the singularities of the linear ODE of x*-'^\ are obtained from this 
modular Atkin- polynomial '89] (j56p together with the covering (j43p (or (|52p ). 

These results strongly suggest that a// </ie singularities, and the associated 
polynomials with integer coefficients we obtained in previous papers (from some 
involved Landau singularities analysis [71|, I60j). should have an interpretation as 
singularities of an (absolute) Klein modular invariant j{N ■ t), or equivalently of 
the modular Atkin-polynomials [89j . for higher values of iV, when rewritten in the w 
variable. 

The diflicult part, here, corresponds to find the well-suited covering (|43|) . or ([52)1 . 
"wrapping" the seven singularities in Z2 onto the three singularities 0, 1, 00 of the 
hypergeometric function 2F1. Noticeably, the well-suited covering (|43|) (or ([52]) ) does 
not correspond to a partition according to the critical exponents (|24|) (see section 
(|4.ip '). but to the following selected partition: 

w = 0, w = CO 0, 

- = -lA - 1, (57) 

vj — 1, w — 1/4, 00. 

To find [SO] this covering ^ among the J^l^^ (^) • (2^ - 2) = 3^ - 3 • (2^ - 2) - 3 = 
1806 possible ones, one had to see the seven singularities on the same footing. 

J In the case of modular Atkin-polynomial of degree one in j, like II56I I. one sees immediately the 
rational parametrization l|51|l . The introduction of the modular Atkin-polynomial becomes necessary 
when one does not have a rational parametrization like (1511 anymore, because the genus of the 
modular curves is no longer zero. 
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6. From 2F1 with a pull-back to 3F2 with a pull-back. 
6.1. Linear differential equation for Xd\t). 

For Xd \ (see (HZD), we choose x = t ' as our independent variable. We find that 
the linear differential operator for {^) is of order six, and has the direct sum 
decomposition 

^ ©4') ©43), (58) 



with 



if) = Dx + -^Hx-l), 
dx 

Lf' - Dx- ^2 , . Dx + , (59) 



L^^' = Dx 



{l + x){x-l) {l + x){x-l)x' 

(3) _ 

3 (8 + 36 + 63 a;"* + 62 + 21 a;^ - 6 a; - 4) 
2 (a; + 2) (1 + 2 a;) (1 + a;) (a; - 1) (1 + a; + x"^) x 

ni{x) 

(a; + 2) (1 + 2 a;) (1 + xf (1 - a;)^ (1 + a; + x2) a;^ 

n^jx) 

(a; + 2) (1 + 2 a;) (1 - a;)3 (1 + a; + x^) (1 + x)^ x"^ ' 



Dx- 



Dx 
(60) 



where: 



no(a;) = 2 a;^ + 8 - 7 .t^ - 13 a;^ - 58 x'^ - 88 a;^ - 52 - 13 a: + 5, 
nx{x) = 14a;® + 71a;^ + 146a;S + 170a;^ + 38a;^ 
- 112^3 - 94a:^ - 19a; + 2. 

(3) 

The linear differential operator of order two, is equivalent to the second order 
operator corresponding to the complete elliptic integral E(x): 

Ce= Dl +^ + ^7—^. (61) 

X (1 — a;) (1 + a;) 

(3) 

Consequently, this order- two linear differential operator L2 is globally nilpotent: 
actually we have calculated its p-curvature and found that the corresponding 
characteristic polynomial (or minimal polynomial) reads 

(3) 

Similarly, the order-three linear differential operator L3 ' is globally nilpotent: 
we have calculated its p-curvature and found that the corresponding characteristic 
polynomial (or minimal polynomial) reads T^. 

6.2. The solution of as a 3F2 with a Hauptmodul pull-back. 

(3) 

The order-three linear differential operator corresponds to a generalization of the 
weight- 1 modular form operators which have been discovered in the previous section 
([S]). Actually, introducing the rational function 

_ {l+2x){x + 2) 
' ~ {l-x){l + x,+x-y ^^^^ 
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one can find that the three solutions of Lg ' are two MeijerG functions 

p{x) ■ MeijerG ([[-1/2, 1/3, 2/3], []], [[0,0], [0]], Q) , (63) 
p{x) ■ MeijerG ([[-1/2, 1/3, 2/3], []], [[0, 0, 0], []], -Q) , 

and a 3^2 hypergeometric functiorQ: 

p{x)- 3F2 ([1/3, 2/3, 3/2], [1,1]; g) (64) 

where Q is nothing but the Hauptmodul (the reciprocal of the Klein's invariant) with 
the (elliptic lambda function) x changed into — x: 

27 [1 + xfx^ ^ 1 ^ , . 

^ 4 (x2+x + l)3 J{~x) 3{~x)' ^ ^ 

With ([64|) we have a straight generalization of the previous weight- 1 modular 
form where the expression for he was given as a 2F1 hypergeometric function with a 
pull-back corresponding to the Hauptmodul ([5^ . One has here a 3F2 hypergeometric 
functior lffi with a pull-back corresponding to the Hauptmodul ((65)) . If one seeks for a 

modular form interpretation of the solutions of ig^' and, in particular, for (|64p . one 

(3) 

would like to see L3 as equivalent to the symmetric square of a 2^1 operator with 
the same Hauptmodul pull-back ([55]) . In this respect, it is worth recalling the Clausen 
identitj0: 

[a + b+1/2]; zY = (66) 



3F2([2a, a + b, 2 &], [a + 6 + 1/2, 2a +2b]; z^, 

and in particular 

2Fi([l/6, l/3],[l];z)' = (67) 
3i^2([l/3, 2/3, 1/2], [1, l];z), 

which is one of the four classesU of 3i^2 found by Ramanujan [24 that are squares of 
2-Fi representations of complete elliptic integrals. The symmetric square of the second 
order operator 

i^, = 3 



2 x{l + x){x^ +X + 1) '-^ 2 {x^ +X + 1)2 ' 
annihilates the left-hand side of (p7|) where z is taken to be equal to the "Hauptmodul 

^3 



(|65p . We actually found out that this symmetric square is equivalent to with two 



simple order-one intertwinners Vi and V2: 



V2- Sym^^Wi'j ^ if^-Vi, where: (68) 

x-{l + x) 1 2x'^ + 5x + 2 



(l-a;)2 " 2 {x^ +x + l){l-x)' 

t Reminiscent of, for instance, 2i^i([l/3, 2/3], [1]; Q) in (23) of [911. 

ft Note that this hypergeometric function is not a well-poised hypergeometric function. 

^ For such selected arguments od the hypergeometric function we also have the highly remarkable 

Gauss-Kummer quadratic relation: 2Fi{[a, fe], [a-|-6-|-l/2]; 4:z{l-z)) = 2^'i([2a, 2 6] , [a-|-b+l/2]; 2). 

I Four arithmetic triangle subgroups, commensurable with the full modular group, yielding four 

hypergeometric representation of periods of elliptic curves |28| . 
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Equivalently, we could have used the foUowing identitjQ] to rewrite (|64p : 

3F2([l/3,2/3,3/2],[l,l];g) = 2^1 ([1/6, 1/3], [1]; q)' (69) 

+ ^-^Fi ([1/6, 1/3], [1]; q) • 2i^i ([7/6, 4/3], [2]; q) 

where: 2^1 ([7/6, 4/3], [2]; q) = 18 ^^^^i ([1/6, 1/3], [1]; q) . 

This generalization of the Clausen identity concludes our modular form interpreta- 
tion [92j [93] of the third order operator L3 . 

7. Global nilpotence without integral representation: three-choice 
polygons, directed compact percolation, vicious walkers, . . . 

Let us also give here three more examples of global nilpotence that do not correspond 
to n-fold integrals, emerging from the theory of the Ising model, or to n-fold integrals 
of the Ising class, but to enumerative combinatorics. Their global nilpotence, or the 
fact that they could be DFG cannot be immediately seen from a representation as 
an n-fold integral of an algebraic expressior{^. Of course one can always imagine to 
see this global nilpotence as a consequence of a detailed analysis of the corresponding 
series expansions, showing explicitely that they are probably arithmetic Gevrey series 
and in fact G-functions. 

Along this line, it is worth noticing that, for instance, the series expansion of the 
X*-"-* in the variable w are series expansions with mtef^er coefficients (see ((H])). The fact 
that these series of holonomic functions are series with a finite radius of convergence 
and with integer coefficients, gives us a strong prejudice that we are studying G- 
functions. For the various examples of enumerative combinatorics displayed in this 
section, we have a similar strong prejudice in favour of arithmetic Gevrey series: G- 
functions. Instead of such an arithmetic approach, we prefer to consider directly the 
global nilpotence, performing p-curvature calculations modulo large set of primes, 
since this approach is simple, algorithmic and effective. 

7.1. Global nilpotence without (known) integral representation: three-choice polygons 

An order-eight Fuchsian linear differential operator, was found for the perimeter 
generating function of the thrce-choiceil polygons [9^. It is the direct sum of two 
order-one operators and an order-six linear differential operator A/g. 

This order-six linear differential operator Mg, is the product of an order-three, 
an order- two and an order-one linear differential operator, respectively M^, M2 and 
Mi: 

Me ^ M3 ■ M2 ■ Ml, where: Mi = - ln(l - x), 

dx 

t Other rewritings involve the Houn functions //eim(10, 23/9, 1/6, 1/3, 3/2, 1, 1 — Q) and 
Heun{W, -35/18, -1/3, -1/6, 1/2, 1, 1 - Q) • (1 - Q)-^/^. 

^ Recall that the integral of a closed differential n-form over a closed n-cycle is said to come from 
algebraic geometry. In our case we have an algebraic n-form, so we automatically have a closed 
n-form and our integration on hypercubes and not cycles, can be reduced to cycles because the global 
nilpotence property |94) remains stable |95l by the required extensions (desingularization, relative de 
Rham cohomology, etc.)[Y. Andre e-mail]. 

§ Note that there is a factor 4 in the definition of x between the ODE displayed here and the one 
given in | 96| . 
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^2 I -^11 7-. I ^21 



^'^2 = Di + -ii- D., + (70) 

(x - 1) • Pi2 • a; 4 • P22 

where the two order-three and order-two operators are given in (jB.ip and have, beyond 
apparent singularities, singularities at cc = and x = 1 and at the roots of the 
quadratic polynomial + Ax -\-7 . Let us consider the order- two linear differential 
operator M2. The indicial polynomials of these singularities read: 

a; = 0, r-{r + 2), 1% + Ax + 7 x'^ = 0, r-(r-l), 

X = I, (2r-f3)^ X = 00, (2r -f 1) • (2r - 1). 

All the formal series around these singularities have logarithms, except a; = 00 which 
has only square roots singularities {x~^^'^). Keeping in mind that the wronskian of 
M2 is a rational function, one deduces [84] that the differential Galois group of M2 is 
S'L(2, C). 

We have calculated the p-curvatures of the two linear differential operators of 
order three and two of (jB.ip . We found that M3 (resp. M2) is globally nilpotent, 
the characteristic polynomial of its p-curvature being (resp. T'^) for almost all 
primes. 

7.2. Global nilpotence without (known) integral representation: directed compact 
percolation 

Generating functions, associated with the mean cluster size and length for the directed 
compact percolation problem, satisfy various linear ODEs that are displayed in |97) . 
For instance, equation (12) of [97], when rewritten in a homogeneous way, corresponds 
to the order-three linear differential operator which is the direct sum Li © L2, of an 
order-one operator Li and an order- two operator L2" 



dt V2i4 + 2i3-6i2 + 4i-l 
2 12^* - 12<'^ - lli^ + 8i- 1 

- ^* {l-2t) {t-l)t{l + 4t-At^)'^' 
Af^ -12t^ + 7t^ -2t + l 
{l-2t)^l + At-4t^){t-l)^t^ ^ ' 

Similarly equations (16), (18) and (21) of ^97j, when rewritten in a homogeneous 
way correspond to the order-three linear differential operator which is the direct sum 
© of order-one operator L'^ and an order- two operator L'2. For instance for 
equation (16) of [57] we have: 

^1 = +dtHAt^+8t-l. 



24^^ - 72t^ + 46t^ + 8t - 7 
2t) {t-l) (1 + 
u4-24t3 + 2<2 + 28<-9 



^2 ^* {l-2t){t-l) {I +At -Af^)t'^' 



' (i-l)(l-2t)(l +4< -4i2)t2 C^^) 

All these order-two operators (L2, L'2, ...) can be seen to be equivalent and are 
globally nilpotent. 
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7.3. Global nilpotence without (known) integral representation: vicious walkers 

Another example of enumerative combinatorics corresponds to the vicious walkers and 
friendly walkers generating functions [HE]- Equation (4.34) of [HH], when rewritten in 
a homogeneous way corresponds to the order-three linear differential operator which 
is the direct sum Li ® L2, of an order-one operator Li and an order-two operator L2: 

d / t^ 



dt V 1 - t + 3 12 
. 16i2 + 21t-4 4^2 + lot- 3 

Equation (4.38) of [98j . when rewritten in a homogeneous way, corresponds to the 
order three linear differential operator which is the direct sum L[ (S) L'2, of an order- 
one operator L'l and an order- two operator 



dt V(t- 1) (1 - 2<-|-4i2) 



, 5-35t+16t^ 

L'., = Df +2 — — • Dt 

^ * (8t- 1) (t - l)t 

5- 18t-41t2 + 10t3-(-6i4_4t5 
^ (l-8i)(l-02(l + i)'i' 

where these last two operators (L2 and L2) of order two are equivalent and actually 
correspond to a Heun functior[§|. A solution of Li Q) L2 corresponds to a generating 
function with integer coefficients namely: 

^ • (i7eun(-l/8, -1/4; -1, -2, 2, -2; -t) -{l-t + 3t^] 

= 1 +2t +6t^ +22t^ +92t'^ +'i22t^ + ■■■ (75) 

We have calculated the p-curvature of the corresponding linear differential operators 
and found that they are all globally nilpotent. 

A large set of other enumerative combinatorics problems for which a linear ODE 
exists for the generating function, but a representation of that generating function as 
a multiple integral has not yet been found, could be listed [551 llOOl 11011 11021 1103j . 



8. Beyond order two and three ODEs: staircase polygons and Calabi-Yau 
type ODEs. 

Krammer's counterexample to the first Dwork's conjecture (reduction of global 
nilpotence to 2^1 hypergeometric function up to pull-back by a rational function) 
corresponded to a selected Heun function. Let us recall that some 3F2 hypergeometric 
functions for selected values of the parameters [104j like 

3F2([a, /3, 7 + l],['5, 7];a:) (76) 

§ Note that the notation in |98l llOOl are Snow's notations which corresponds to a change of sign of 
the second argument in the Heun functions compared to maple's notations. Also note a sign misprint 
in 1981 {t in 1981 instead of —t in llOOl ). Do note that changing —1/4 into -1-1/4 in the previous 
Heun solution, one gets Heun{— 1/8, 4-1/4, —1, —2, 2, —2; —t)/t^ which is solution of a second order 
operator almost identical to L2 in II73I I. where 4t^ + lOt — 3 is replaced by 4 1-^ + lit — 3. This 
new operator is not globally nilpotent, it differs from L2 by the so-called [501 1511 1521 153) "accessory 
parameters" . 
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have a Heun representation. We want, now, to explore globally nilpotent higher order 
linear operators beyond the 2-^11 Heun and framework. 

In [Appendix D[ the linear differential operator for Xd"^ (0 is displayed and 
factorized in direct sums of an operator of order one, an operator of order three 
and an operator of order four. The operator of order three is seen to be equivalent 
to a symmetric square of a second order operator corresponding to a hypergeometric 
function 2-F'i- Similarly to the situation described in section (|6.2p . the order- four 
linear differential operator is not obviously reducible to a symmetric cube of a 2F1 
(even with a pull-back), it may however be equivalent to a symmetric cube of a 2F1. 
It may well be a 4^3 hypergeometric function with an involved (Hauptmodul, . . . ) 
pull-back, but, for the moment, we have not been able to decipher this order-four 
operator. Let us display, in the following, a few more order-four (and more) globally 
nilpotent operators, related to staircase polygons and Calabi-Yau manifolds. 



8.1. Staircase polygons. 

In an enumerative combinatorics framework other Heun functions can be found for 
the generating function of the staircase polygons [87l I105j . 

The Fuchsian differential equations corresponding to the staircase polygons 
generating functions in d dimensions, that we denote Zd are given in [57] : 



^ a; • (1 -x) (1 - 9 a;) ' " x ■ {1 - x) {1 - 9 x)'' 

7 _ n3 o (l-30a;+128a;^) . 
- +^x.(l-4.x)(l-16x)'^- 

(1 -68a; + 448a;2) 1 



(l-20x-f 27a;^) l-3x 



x2-(l-4.T)(l-16a;) x^-{l-4:x) 

and we give Z5 and in [Appendix B.5[ 

All these linear operators (up to Z^ in [87]) are globally nilpotent: this is a simple 
consequence of the fact that the generating functions of the staircase polygons are 
expressed as n-fold integrals of (very) simple algebraic integrands, and are therefore 
DFG. We calculated systematically the p-curvatures of all these Fuchsian linear 
differential operators (considered modulo the first thousand primes) and found that 
their characteristic polynomial is T'^~^, where d — 1 is the order of the operator Z^. 

The first two staircase linear differential operators Z3 and Z4 actually correspond 
to Heun functions. The solution of the Heun linear differential operator Z3 and its 
series expansion has integer coefficients: 

iJeun(l/9, 1/3,1,1,1,1,.t) = i/eun(9, 3, 1, 1, 1, 1, 9 x) = 

= l + 3x+l5x^ + 93x^ + 639x* + A653x^ -f 35169x^ 

The iJeun(l/9, 1/3, 1, 1, 1, 1, x) solution can also be written in terms of a modular 
hypergeometric function (corresponding to the weight- 1 modular form Iiq in [86] ): 

iJeun(l/9, l/3,l,l,l,l,a;) = (77) 

L/ 

{x-l){9x-lf 



(3 a;- 1) • (243 a;^ - 243 + 9 a; - 1)) • 2i^i( [1/12, 5/12], [l],Me) 



with: Mfi = 1728 



(3 a; - 1)3 • (243 a;^ - 243 a;2 + 9 a; - 1)3 
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Aa ^ Di +( - + + — • D 



Joyce has shown ^79 that the square of this Heun functioi|j| is related (up to 
quite involved algebraic transformations of the arguments see also equations (23), 
(24) in [57]) to the simple-cubic Green function where its DFG nature becomes clear: 

J J Jo 1 - z/3 • (cos(a;i) + cos(a:2) + cos(a:3)) ■ 

This Heun function Heun{l/9, 1/3, 1, 1, 1, 1, x), we now see as a weight-1 modular 
form, had been seen by Guttmann and Prellberg to be equal to the product of 
two complete elliptic integrals of the first kind K (up to slightly involved algebraic 
transformations namely equation!^ (25) and (28) in [87|). As a (curious) byproduct 
one thus deduces a new relation between a weight-1 modular form and a product of 
two complete elliptic integrals. 

The function Heun{\/A, 1/8, 1/2, 1/2, 1, 1/2, Ax) corresponds to the Heun linear 
differential operator: 

2-(8a:-l) 

^x ' 4a;-l ' l&x-lJ ' (4a;- 1) (16x- l)a;' 

Its series expansion has integer coefficients: 

iJeMn(l/4, 1/8,1/2,1/2,1,1/2, 4a:) = 

= iJeun(4, 1/2,1/2,1/2,1,1/2, 16 a;) = 

= 1 + 2 X + 12 a;2 + 104 .T^ + 1078 x^ + 12348 a;^ -I- 150528 x^ + 

One easily verifies that the symmetric square of the previous Heun linear differential 
operator is nothing hut Z^. Actually 7?eMn(l/4, 1/8, 1/2, 1/2, 1, 1/2, 4x)^ is 
solution of Z4. 

To our knowledge this Heun function iJeuri(l/4, 1/8, 1/2, 1/2, 1, 1/2, 4a;) does 
not have an interpretation as a weight-1 modular form (see also [SHI 11061 11071 [55]). 
With these two Heun functions we are thus still in a "Dworkian" framework 
(hypergeometric functions up to a pull-back, selected Heun functions, etc.). 

The very nature of Z5, Z^,, Z-j, . . . , is less clear. In particular, it is far from clear 
that they can be written in terms of 2F1 hypergeometric functions. Let us introduce 
some kind of "multi-singular" (beyond the four singularities of a Heun ODE) second 
order linear ODE: 

^^~^^\x 3(1 -a;) l-9a; 3(1 - 25a;) J " 

19 171 5875 49 



270(l-a;) 40(1 -9x) 216(l-25a;) 30a; 
1 9 125 

^ 90(l-a;)2 + 10(1 -9a;)2 ^ 18(1 -25a;)2 ^ > 

One finds that the explicit expression of Z5 is "close" to be a symmetric cube of 
this "multi-singular" (beyond the four singularities of a Heun ODE) linear differential 
operator: 



Z5 = Sym (zs) + — 



100 (1 - xY (1-9 xY (1 - 25 xY ■ x3 ' 

X See also equation (36) in 98 for the generating function of a watermelon counting (union of friendly 
walkers) . 

§ Remarkably, this is also true for _fi'eMn(l/4, 1/8, 1/2, 1/2, 1, 1/2, 4t)^ (see below) (these are 
equations (25) and (28) in (STj). 
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where P{x) is a polynomial of degree eleven with integer numbers. In fact this order 
four operator is associatecQ with Calabi-Yau manifolds [108j (see also next section). 

Similarly, introducing the "multi-singular" (beyond the four singularities of a 
Heun ODE) second order operator 



V ,1 1 4 9 , 
ze = Di+(-- — - — - — • 



1 - Ax 1 - 16a; 1 - 36a; 
1 p{x) 



(79) 



5 a;-(l-4a;)2(l-16a;)^(l-36a;)^' 
p(x) = -7 + 696 a; - 22224 a;2 + 298816 a;^ 
- 1603584 a;"* + 3068928 a;^ 

one finds that Zg is (up to some rational functions A{x) and B{x)) "close" to be a 
symmetric fourth power of this previous "multi-singular" (beyond the four singularities 
of a Heun ODE) operator ((79|) : 

Ze = Sym^izfi) + A{x) ■ + B{x). (80) 

Let us recall that, even for a second order operator, it is not so easy, and 
systematic, to see that a solution reduces to a hypergeometric function 2F1 up to a 
possibly involved rational pull-back and up to a multiplication by some iV-th root of 
a rational function (see ((37)) . (|49|) ). It has to be the case [83] for zero-curvature second 
order operators (and this is an important motivation to perform our p-curvature 
calculations, to extract, very quickly, the zero-curvature situations like (I37|) ). when 
we have a basis of algebraic solutions, but this is far from clear, at first sight, for a 
globally nilpotent operator (see (P5|) ). 

When one encounters a higher order irreducible linear differential operators like 
these Zrf for d > 5, (or other order-three irreducible operators like M3 in the three- 
choice polygon problem, see (jB.2p ). it remains to see|3| if these irreducible operators 
cannot be symmetric powers of a smaller order operator with an emphasis on second 
order operators simply reducing, or not simply reducing, to 2F1 hypergeometric 
functions. 

With these globally nilpotent staircase operators Z5, Zq, Z7, . . . we are clearly 
in a win-win situation. In a first scenario, the very simple and regular accumulation 
of singularities (a; = 1/'^^, for n even or odd) cannot be wrapped onto three, or four, 
canonical singularities 0, 1, a, 00, and one gets brand new examples of selecteci§| 
"multi-singular" (beyond the four singularities of a Heun ODE) globally nilpotent 
operators of arbitrary order d—1, and fascinating algebraic geometry interpretations 
(similar to the one Krammer gave for the Heun ODE like ((S]), with Shimura curves) 
which remain to be found. In a second scenario these staircase generating functions 
can be expressed as products of complete elliptic integrals (up to involved algebraic 
transformations) or 2^1, 3^2, with some (Hauptmodul . . . ) puUback, and this would 
provide a set of highly non-trivial effective algebraic geometry results associated with 
these n-fold integrals and these enumeration of staircase polygons. 

f A. J. Guttmann, private communication. 

^ To check if an order q operator is exactly a symmetric power of, say, an order-two operator can be 
done systematically (pattern matching). To see if an order q operator is equivalent (homomorphic) 
to some symmetric power, is less easy and systematic 11091 IllOl 111 II 11121 1113| . 
§ Do recall that, for instance, Heun operators are not generically globally nilpotent (see l|8)l ). 
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8.2. Calabi-Yau type ODEs. 

Other non-trivial examples of non-trivial globally nilpotent high order operators 
can be found with the fourth order differential equations of the so-called Calabi- 
Yau type |114[ I115j . We will not give any detail on the construction of these new 
type of Calabi-Yau manifolds using conifold transformations from toric Calabi-Yau 
hyper surf aces. We just display some of these fourth order differential operators. 

A first order-four differential operator comes from Kontsevich's observation |114j 
that two selected matrices for a quintic and its mirror, actually correspond to 
monodromy matrices of the Picard-Fuchs operator: 

~ 5'x-{e +ljie +^){e + ^)ie + ^), (si) 

where: 9 ~ x ■ — 

ax 

having four solutions which can be expressed in terms of hypergeometric functions 

/ n n n n n + ii n + 22 n + 1 \ 

'Hh' 5' 5' 5^J' 
where n — 1, 2, 3, 4 and «i, Z2, «3 are three integers in the list [1, 2, 3, 4]. Another 
example of order-four differential operator which can be found in 114j reads: 

- X- {659^ + 1309^ + 10502 -|- 400 + 6) 

+ Ax^ {4:9 + 3)i9 +lf {4:9 + 5). (83) 

The critical exponents of ((83|) are given in [ll4j in P-Riemann function notations. 
They are very simple integer or rational numbers (3/4, 5/4). 

Associated with the diffeomorphisms [115) Xfl^ ^20 another order-four differential 
operator reads: 

9'^ - 2 a; (102 6** -1- 204 6'^ + 155 9'^ -\- 53 9 +7) 

-I- 4x2(61 + 1)2 . (396 + 792 61 +311) (84) 

- 784x^ {9 + l)(6i + 2)(26i + l)(26i + 5), 

which, simply written in x, has a form very similar to the staircase operators Z5 given 
by (|B.4[) in [Appendix B.5[ namely : 

4 1568x2- 268^+1 3 

" x(l-4a;)(l-196a;) " 

7 - 2962 a; + 39260x2- 116816x3 , 

D 



(1- 196x)(l-4x)2x2 
1 - 1028 X + 22740 x2 - 90944 x^ 

(l-196x)(l-4x)2x3 
„ 7 + 3920x2 - 622 X 



(85) 



x3 • (1- 196x)(l-4x)2' 

Associated with the diffeomorphisms 115^ ^200 140 another order- four linear 
differential operator reads: 

9* - X (113 61'* + 226 9^ + 173 6*2 + 6O6I +8) 

- 8x2 (gi^_i)2 . (^iigQ2 2389 + 92) 

- 484 x^ {9 + l){9 + 2){2 9 + l){2 9 + 5), 
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which, simply written in x, has, again, a form very similar to the staircase operator 
2^ given by (|B.4p in [Appendix B.5[ or to the previous order four operator (|85p . but 
with the denominators (1 - 196 x) (1 - 4 a;)'' of ((511) changed into (1-121 (1 + 4 xy. 

We have calculated the p-curvature of these four order-four linear differential 
operators ([STjl , ([83| , (|85p and ([86|) (modulo the first thousand primes) and found that 
the characteristic polynomial, as well as the minimal polynomial, read T*, and that 
the Jordan block reduction of the 4x4 p-curvature matrix reads: 

"0100' 
10 

(86) 

1 




9. Beyond holonomic functions: ratio of holonomic functions. 

In [1161 I117j it has been shown that the enumeration of three-dimensional convex 
polygons can be written as ratio of holonomic functions. 

For instance the equation in Proposition 4.12 in [117 , or the equation just before 
the conclusion in 116|, gave the perimeter generating function for three-dimensional 
oriented convex polygons as the ratio 

— i-, with: Ni= Ai{u)-S3{x)+ A2{u), (87) 

^3(X) 

where u denotes the square root of x, where Ai{x) and A2{x) denote algebraic 
expressions of u and where S3{x) denotes the solution of the previous staircase 
operator Z3 given by (j77p . Such a numerator Ni, is also of the type considered 
in section (O, namely functions for which an integral representation has not yet been 
found. Do note that ratio of holonomic functions are far from being holonomicjj, as 
can be seen on the solutions of the Chazy III non-linear ODE (see for instance (3.37) 
in |118| . and page 1878 in 4J). In terms of x and not its square root u, the algebraic 
expression A2 is solution of a linear differential operator of order four, L4 (direct sum 
of an operator of order three and an operator of order one). The algebraic expression 
Ai is solution of a linear differential operator of order five, L5 (direct sum of an 
operator of order two and three operators of order one), the product Ai(u) ■ S^lx) 
being solution of a linear differential operator of order ten Liq (symmetric product 
of and is). The numerator A^i is solution of an order fourteen linear operator 
that can be obtained as the LCLM of Liq and L4. We have calculatec0 the p- 
curvature for the first two hundred primes, of this order fourteen operator (the degree 
of the polynomial coefficients is 320) and found that its characteristic polynomial is 
T^'^, its minimal polynomial being T^. The global nilpotence of this order-fourteen 
operator is a straight consequence of the expression of A^i in (jST]) . Examples like 
the enumeration of three-dimensional convex polygon suggest to seek for new classes 
of solutions that arc, not only ratio of holonomic functions, but ratio of solutions of 
two globally nilpotent operators "algebraically equivalent" in the sense of relation (|87p . 

tt In contrast the product of two holonomic functions is holonomic. 

^ Note that these p-curvature calculations are very quickly performed, when the factorization (or 
LCLM-factorization) of such order fourteen operators is easily reaching the limits of our 32-Gigas 
computer facilities. 
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such ratio of periods probably having a (modular) interpretation as r functions and 
Painleve-like Picard-Fuchs deformations [119j . 

10. Beyond global nilpotence: linear differential operators with irregular 
singularities 

We have encountered in previous publications ^75j n-fold integrals annihilated by 
(minimal) linear differential operators which are obviously globally nilpotent, namely 
the two-point correlation functions and the form factors of the off-critical lattice Ising 
model. The linear differential operators of the form factors have a nice "russian-doU" 
structure (see ((T8| V The linear differential operators Fj{N) occurring [75] in these 
factorised russian-doU form were seen to be equivalent to symmetric powers of the 
second order linear differential Le corresponding to the complete elliptic integral of 
the second kind E. Consequently, they are obviously globally nipotent. The scaling 
limit of these linear differential operators also exhibit a "russian-doU" structure, but 
they are not globally nilpotent. 

(n) 

The scaling limit of the fj^ ^'s amounts, on the functions, and on the 
corresponding differential operators, to taking the limit — > oo and t — > 1, 
keeping the limit x = N ■ {1 — t) finite, or in other words, to performing the change 
of variables t = 1 — x/N, keeping only the leading term in N. Performing these 

(n) 

straightforward calculations, the linear differential operators in t for the fj^ j^'s where 
N was a parameter, become linear differential operators in the scaling variable x. 
Calling the scaling limit of the operator Fj{N) we found [75] for j odd, 

that 

ipscal T seal 

^1 = ^2 ^ 

rpscal r seal j seal 

^3 — ^4 ■ ^2 ^ 

Fiscal T- seal j seal r seal 

5 — ^6 ■ ^4 ' ^2 ^ 

where 

Lf^' = f 6 Dl + 96 x^ Dl + A0{2-x^) x^ DI 

+ 8{x^ ~2)xD^ -I- 9 a;"* - 8 -t- f 6, 

L^"''' = 4x-Dl +AD^ -X, (88) 

and LIq°'\ Lg™', Lg^"' are given in ^75j. Similar relations occur for j even 'TP. Thus, 
we see that the scaled operators F^'^"-'' have a "russian doll" structure inherited from 
the lattice operators Fj{N). 

Consider the linear differential operator corresponding to the modified Bessel 
function Bessel{n, x/2) for n — 0, namely: 

B ^ Dl +^ (89) 

We recognize, in this linear differential operator, the exact identification with the 
scaled differential operator = L?^"''' . We find that the symmetric square of 

the linear differential operator B, and the scaled operator L^"''' are equivalent, the 
symmetric third power of the linear differential operator B, and the scaled operator 
j^scai g^j.g equivalent, and, more generally, the symmetric j-th power of (|89)) and the 
scaled operator L*™/ are equivalent, Lj'^j ~ Sym-'(i3). 

Global nilpotence implies fuchsianity. The scaling limit generates a confluence of 
the regular singular points [1201 1121| we had on the lattice, yielding linear differential 
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operators, which are not Fuchsian anymore because of an irregular singular point at 
infinity: we are leaving the universe of G-functions for the universe of "Hamburger" 
functions [21 . 

Let us explore, however, the p-curvatures of the previous non-Fuchsian hnear 
differential operators which correspond to scaling limits of globally nilpotent linear 
differential operators. The calculations give, modulo the prime p, the following 
characteristic polynomial for Lj^"': 

and the following characteristic polynomial for Ll"^"': 

T+P—^) . fr+^) . + ■ (t+P-±^). (91) 



We have also calculated the p-curvatures of and the corresponding characteristic 
polynomials. For almost all primes, this characteristic polynomial has a very simple 
expression one can see as a simple deformation of the characteristic polynomial of a 
nilpotent operator: — T — (T — 1) • T • (T + 1). Similar calculations performed 
for Lg'^"', Lg'^°', L^'^"', up to Lf™' give the following results for the characteristic 
polynomial of the corresponding p-curvature: 

T'^^) ■■■ {t-2] ■■■ (t-i) X (92) 

n-1 



2 

T- (T+1) ••• (T + 2) ••• (T 



for Lf^'^' with n odd and 

i = n/2 



n 



T+^+z], (93) 



1=1 -n/2 

for L^^°' with n even. All these calculations have been performed for all the primes 
p < 100. 

A remarkable structure "beyond global nilpotence" clearly remains to be 
discovered by mathematicians for the functions of "HamburgeiQ type" (one irregular 
singular point at infinity) that occur in field theory, or, more simply, in the scaling 
limit of DFG holonomic functions of lattice problems. 



11. Conclusion 



One can probably conjecture that when the generating functions of the various 
problems of enumerative combinatorics are found to be solutions of Fuchsian ODEs, 
quite systematically the corresponding linear differential operators are globally 
nilpotent, these holonomic functions being "DFG", their rewriting in terms of n- 
fold integrals being just a question of time, work and/or stamina. The generating 
function of the perimeter three-choice polygon, of the directed compact percolation or 
of the vicious walkers are such examples. In this paper, we have studied a quite 
large number of n-fold integrals of algebraic integrands and their corresponding 

X To sum up quite brutally the situation, one may say 1191 1201 that almost all the special 
functions occurring in theoretical physics are either G-functions for Fuchsian ODEs or "Hamburger" 
functions 1211 when an irregular singularity occurs. 
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Fuchsian ODEs. In particular, we looked at the p-curvatures of their factors, not 
to see if these linear differential factors were globally nilpotenfl, but to understand 
how these differential factors "succeed" to be globally nilpotent. One must keep in 
mind that the Fuchsian ODEs for the n-fold integrals of theoretical physics (Feynman 
diagrams 63j , • ■ ■ ) are generically of quite large orders (as an example the minimal 
order ODE for x*^^^ is of order 33, see [5]). Since the corresponding minimal order 
ODEs are necessarily globally nilpotent because they are DFG, the question one can 
ask is how an order 23, 33, 50 linear differential operator succeeds to be globally 
nilpotent? Throughout all the examples displayed in this paper we have seen that 
the (minimal order) linear differential operators of quite large order actually factorise 
into products (sometimes direct sums and products) of linear differential operators 
of smaller orders (one, two, sometimes three and four). The global nilpotence of the 
order-one operators just corresponds to wronskians that are A^-th roots of rational 
functions, most of the order-two linear differential operators and, for instance, the 
order-three and four operators occurring in the factorizations of x^^^ and x^^^ j having 
a typically "Dworkian" interpretation since they correspond to either the second order 
linear differential operator associated with the complete elliptic integrals (of the first 
or second kind), or equivalently to 2F1 hypergeometric function or to symmetric 
powers (square and cube) of these second order linear differential operators. More 
remarkably other second order linear differential operators were found to correspond 
not only to globally nilpotent operators but to zero curvature operators. The solutions 
are algebraic functions corresponding to selected algebraic curves. For instance, we 
encountered genus six and genus five algebraic curves (for and <I>^^ respectively) , 
their roots being expressed as complete elliptic integrals of the third kind with a 
"characteristic" corresponding to genus three curves. As far as algebraic curves are 
concerned, we also gave many examples of zero curvature linear differential operators 
of different orders corresponding exactly to algebraic functions associated to modular 
curves, namely the A-extensions C{N, N; A) for some selected values j75j of A. From 
a "Dworkian viewpoint" let us recall that these algebraic functions can be expressed 
as pull-backs of 2F1 hypergeometric functions. We also gave many examples of 
zero curvature linear differential operators of different orders corresponding exactly 
to algebraic functions associated to modular curves. As far as second order linear 
differential operators are concerned, the most spectacular example came from the 
linear differential operator Z2 occurring in the factorization of the linear differential 
operator for x^^''- The global nilpotence of that operator Z2 was seen to correspond 
to a highly non-trivial pull-back of 2F1 hypergeometric function, namely the weight- 1 
modular forms /ig ! Most of these examples correspond to n-fold integrals associated 
with the Ising models or more generally n-fold integrals of the so-called Ising class. 

Among all the globally nilpotent operators we displayed in this paper, other n- 
fold examples came from enumerative combinatorics and others from Picard-Fuchs 
(Gauss-Manin connection [7]) constructions on (mirror) Calabi-Yau hypersurfaces 
with conifold singularities. In the first case (enumerative combinatorics) only two 
examples were clearly "Dworkian" the n-fold integral being expressed alternatively as 
a Heun function or as a pull-back of a 2F1 hypergeometric function. All the other 
examples seem to go beyond a strict hypergeometric "Dworkian" framework. One 
seems to explore, similarly to Krammer counterexamples, Heun functions that cannot 

K This question has been solved by mathematicians: these n-fold integrals are holonomic functions 
with rational critical exponents, and are even DFG. 
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be reduced to 2^1 hypergeometric functions, and, more generally, holonomic functions 
corresponding to linear differential operators of order two, three, four, . . . , with, at 
first sight, many more singularities than 0, 1, and 00. 

The details and the richness of the situations, and structures, encountered with 
our physical examples were certainly not obviously expected from the DFG diagnostic. 

Understanding how linear differential operators are globally nilpotent led us to 
discover different structures on various algebraic varieties (elliptic curves and complete 
elliptic integrals, algebraic curves that are modular curves, weight-1 modular forms, 
. . .) that provide a deeper understanding of the underlying mathematical structures 
"hidden" in the physics problems we study. When we see that, to sum up things 
brutally, the global nilpotence of the linear differential operator for x'"^^ is inherited 
from the global nilpotence of Z2 which corresponds to the weight-1 modular forms hg, 
one understands the "complexity" of the holonomic function x^^\ totally, and utterly, 
differently. An interesting generalization of the previous weight-1 modular form was 
found with and equation (|64p . with the occurence of a 3F2 hypergeometric 
function with a Hauptmodul pull-back. This reinforces the viewpoint we tried to 
promote that a (serious) theory of the Ising model embed all the theory of elliptic 
curves (modular curves, modular forms, . . . ). The last Calabi-Yau examples confirm 
that viewpoint: discovering the underlying algebraic varieties (or projective spaces 
minus some singular sets [1221 1123) ). is fundamental. What are these algebraic 
varieties curves, surfaces, and higher dimensional algebraic varieties? Going a bit 
further beyond the simple mantra "it is derived from geometry" , we tried, and often 
succeeded, to find explicitly the structures of effective algebraic geometry that are the 
"deus ex machina" of our theoretical physics problems. 

One can now propose the following systematic, and quite algorithmic, study 
of every n-fold integral of algebraic expression encountered in theoretical physics: 
first generate large series expansions of these n-fold integrals to find out the 
linear differential operators that annihilate these series, then get the minimal order 
differential operators, then factorize and LCLM-factorise, as much as possible, these 
minimal order differential operators, then examine the irreducible factors to see if they 
are not equivalent to symmetric powers of smaller order operators, then calculate the 
corresponding p-curvature^^ of all these smaller order irreducible operators to see if 
they have zero curvature, or if they are globally nilpotent, and, finally, examine all 
these smaller order irreducible operators to find out if they correspond to 2F1, 3F2, 
4F3 , . . . hypergeometric functions up to a rational pull-back. 

Saying that theoretical physics should eventually reduce to classification of 
singular varieties [122|, I123j , thus reducing most of it to effective algebraic geometry, 
is certainly a too drastic simplification. With section (fTO|) we see that some 
nice generalisation of the notion of global nilpotence does exist, and needs to be 
explored, for the Hamburger functions |21| that typically occur in physics (linear 
ODEs with one irregular singularity, for instance at 00). Let us just say that an 
effective algebraic geometry viewpoint of lattice statistical mechanics, enumerative 
combinatorics, particle physics, solid state physics, theoretical physics, hopefully 
yielding the emergence of a new "Algebraic Statistical Mechanics" is certainly a step 
in the good direction. 

t|For operators of quite large order and degree of the polynomial coefficients, the p-curvature 
calculations are much quicker than the factorization and LCLM-factorisation. 
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Appendix A. Factorisations of multiple integrals linked to ({3) 

In Apery's proof of the irrationality of C(3) a crucial role is played by the linear 
differential operator |124| : 

(t^ -Mt +l)-t^-Df + (6<2 -I53i + 3)-t-D^ (A.l) 

+ {7t'^ ~112t +!)• A +(i-5), 

this operator being linked to the modularity of the algebraic variety: 

1111 

X ^ ^ y + - + z -\ hwH = 0. 

X y z w 

Operator (|A.1|) is, in fact, the symmetric square of the second order operator |125j : 

4- - 34i +1) - i -D^ (A.2) 

+ 4 ■ {2f - bit + 1) ■ Dt + (t ~ 10), 

If one introduces the change of variable 

t = :i:li_M, (A.3) 

1 — X 

the second order operator (jA.2[) becomes 

2 1 - 18 a; n _ 1 10 - llx + 9g;^ 
^+(l-9a;)x " 4 ■ (1 -9x)(l - x)2x' ^ ' 

Gonsidering the solution up to a multiplicative square root (1 — xY^^, amounts to 
perform the symmetric product of this operator (|A.4p with the order-one operator 
Dx + 1/2/ (a; — 1) and transforms (jA.4|l into the second order operator Lx'. 



Dl + — ^^^f--^^ -Dx - 3- 



I -2Qx +21x^ „ , l-3a; 



(1 -9x)(l - x)x (1 -9a;)(l — x) a; ' 

which is nothing but Z3 (see section ()8.ip . or equivalent to the operators Z2 and M2 
of ) that occurred in the staircase polygons, and has the Heun function solution 
(1 — xY^"^ ■ Heun{l/9, 1/3; 1, 1, 1, 1; x) that already occured in Z2 and M2. Again 
we have here [125] the occurence of a Picard-Fuchs equatior0 for the modular family 
of elliptic curves associated to ri(6), that is the weight-1 modular form h^. 

^ Similarly 11261 one also has Gauss-Manin systems of Shimura families of abelian surfaces having 
multiplication by a quaternion algebra over Q. 
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In [5D] we also obtained an order-four Fuchsian linear differential equation (also 
related to the analyzis of C(3)) which factorize^ in four order-one differential operators 
{Dx denotes d/dx): 

2 (3 a; - 1) 



L„ = Dx 



■ Dx" 



+ 



+ 



[x — l)x 

[ix^ +{n'^ +n-b)x - 2n{n + l)) 

(x — 1)2 
{x^ +2n{n + l)) 



■Dx^ 



■ Dx 



[x — 1)2 x'^ 

71 • (n + 1) • ((n^ +n + l)x + (n - 1) (n + 2)) 



(A.5) 



\Dx 



Dx 



{x- 

d\n{Ai) 

dx 
d\n{A3 



Dx 



dln{A2 



dx 



Dx 



dx 
d\n{A4) 
dx 



These order-one linear differential operators have rational solutions. Such factorization 
into order-one linear differential operators with rational solutions is characteristic of 
factorization of globally nilpotent operators when they are considered modulo a prime. 
Here, remarkably, such a factorization takes place for the exact operators over Q and 
not only mod prime ! The Jordan block reduction of the p-curvature of the operator 
(|A.ip read respectively: 

1 








(A.6) 



when it gives matrix Ij86p for (jA.5|) for any integer n. For any integer n the 
characteristic polynomial, as well as the minimal polynomial of the p-curvature of 
(jA.Sp reads T'^. For n a non integer parameter the characteristic polynomial reads 



+ {x-2)P -U -T^ + ((x + l){x- ■ 

'{x 



(A.7) 



with: 



U ^ Ai 



f(x~ 1)2 x^ \ 



An = {n - 1) ■ n ■ {n + 1) ■ {n + 2) ■ ■ ■ (n + p - 2), 
where An, as it should, vanishes modulo the prime p, when n is an integer 



where: 



Appendix B. Display of miscellaneous Fuchsian linear operators of the 
paper. 

Appendix B.l. Operators M2 (resp. M^) for three- choice polygons 

The order-two (resp. three) operator M2 (resp. M3) occurring in the factorisation 
(|70p of the order-six linear differential operator Mg associated with the three-choice 

J Note a misprint in |60| one should read In Ai , instead of Ai , in the equations defining the Ai after 
equation (H.2) in [61] [71] • 
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polygons generating function (see section (|7.ip ) reads 

M2 = Dl + ^ + (B.l) 



with: 



Pii = 3437 a;^- 1341 + 4188 x^- 24160 + 38400 

+ 10752 - 34816 x + 12288, 
P12 = 3437 a;6- 5826 a;5 + 5280 a;"* - 7360x3 

+ 7680 + 3072 X- 4096, 
P21 = 24059 + 35756 x^ + 116792 x^ - 480784 x*' + 693824 x^ 

- 2361856 x^ + 2886656 x^ - 3739648 x^ + 3670016 x - 1376256, 
^^22 = -(16 + 4x + 7x2) . (1 - x)2 . X • p^^^ 



and: 



M3 = P(x) -Dl + ■■■ (B.2) 
where P{x) denotes the (head) polynomiah 

P{x) = 4x3(l + x)(l-x)2(4 + x2) 

X (16 + 4x + 7x2) • p^2(x)3 • g(x), 



Q(x) = 116620x^2 .39739 3,11 ^ 2816770x1° - 4827228 x^ 



- 5350720 x^ + 12343408 x^ + 473056 x^ - 13436096 x 



5 

9007872 x" + 1064960 x^ - 1421312 x^ - 327680 x - 65536. 



Appendix B.2. The order-two operator AI2 in L^o) 

The order-two operator M2 occurring in the factorization L^o) — M3 ■ M2 reads: 
P2- M2 = P2- Dx^ + {x-l)-pi- Dx + 3- Pa, 



with: 



P2 = (x - 4) (x - 2) (x - 1)2 . (2 + x) (4 + 3 X + x^) • x, 
= -64+ 128x+ 196x2 + 20x3-57x'^-14x''^ + 7: 
po = 33x3- 16 - 20x + 44x2-llx'^-9x^+3x^ 



and: 



Ms = q3-Dx^ + ■■■ 

(73 = - 2)2 (x - 1) • x2 • (1 + x)2 ■ p2 -r^, where: 
rg = 1280 - 1344 x - 6848 x2 - 21456 x^ + 82416 x^ + 74876 x^ 



- 44684 x^ - 48873 x^ + 32112 x^ + 25252 x^ + 1728 x^° 

- 1918 x" + 648x^2 _^ 120x^3 -I- 4x^'' - x^^ 

Note that x here is actuaUy Aw. The set of singularities have a w 1/4/^ 
covariance, that is a x <~> 4/x covariance: 

/4\ 4 

P2\ 



(4 \ 4 — X 

x)= ^^^•(T3^-^2(.) 
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Appendix B.3. The order-four operator M4 in L n) 

The order-four operator occurring in the factorization L , (4) = M4 • Ki ■ Z\ reads: 

M4 = q^-Dl+ 93 -Dl +q2-Dl ^qx-Dx ^ Qo, with: 
54 = 16(4-a;)(l-a;)''a;^-Q4, 

Q4 = - 128 - 2233 a; + 2847 a;^- 3143 a;^ + 3601 a;^- 144 a;^ + 64 a;^ 
93 = 32x^1 -xf-Qs, 

Qs = 768 - 4712 + 54621 - 226585 a;'^ + 271255 a;* 

- 253247 x^ + 190228 a;^ - 45848 x - 3328, 
92 = -8a;2(l-a;)2 .Q2, 

Q2 = 23360 x^ - 140752 a;® + 1814065 x^ - 7479930 x^' + 10944040 a;'^ 

- 1 1262286 x^ + 9445431 x^ - 4048776 a;^ + 419280 x + 47104, 
qi = 8x{l- x) ■ Qi, 

Qi = 64640 x^° - 382600 x^ + 5520835 a;® - 22754401 

+ 38212402 x^ - 43444138 x^ + 39867319 x'^ - 22329197 a:^ 
+ 5075028 x^ - 21248 x - 47104, 

go = -65536 + 1444096 x + 4876704 x^ - 79483588 x^ + 250389985 a;^ 
+ 382946518 a;^ - 307163242 a;^ - 380545497 a;^ + 165955737 a;^ 

- 40044089 a;'' + 2440592 a;^° - 419904 a;^^ 

Appendix B.4- The second order operator Ni in L^(5) 

The second order operator A^i occurring in the factorization of L (5) reads: 
Ni = P2-X- Dl + Pi- + x^ ■ Po, 

P2 = -2(l + 8a; + 20a;2 + 15a;3 + 4a;'')(l-a;-3a;2 + 4a;3) 

X (1 + 2 a; - 4 a:^) (1 - 3 a; + a:^) (1 + 2 a;) (a; - 1) (1 + a;) • P2 
P2 = 608 x^^ + 88 x^^ - 3092 a;" - 7329 a;^^ + 2156 a;^^ + 15088 x^^ 
+ 7054 x^° - 3476 x^ - 15856 a;^ - 6198 a;^ + 7668 x^ + 4064 x^ 

- 1014a;*-732a;^ + 36a;^ + 63a; + 8, 

Pi = - 16 -188 a;- 592 a;2- 384 a;3 + 48895500 a:"- 24850920 a;i2 

- 216496824 x^^ + 526808 .t^ - 24630720 x^^ - 203254764 x^^ 

+ 526770352 a;^^ + 243268011 a;^^ + 34759480 x^^ + 91207968 x^^ 
+ 16933696 a;^*' - 13056000 a;^^ - 5926400 a;^^ + 112640 a;^^ 
+ 622592 a;^° + 80791756 a;^^ - 255450647x^2 _ 7415 -j;* 

- 1427696 a;^ - 6112060 x^ + 10047555 .r^" - 686378692 x^^ 
+ 511457760x^0 - 507535276 + 380400540 a;^^ - 33404 a;^ 
+ 113836 

Po = -144- 1344 X- 4416004x^3 - 7073308 x" + 223379312x^2 
+ 107333356 x^^ + 66480 x^ + 1927564 x^ + 2688x2 

- 364322176 x" - 265131272 x^^ + 279826283 x^^ (B.3) 
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- 10176392x24 - 2082592 a;25 + 67584 x^^^ + 311296 a;^^ 

+ 39019332 x^^ + 45440155 x'^^ + 176313 x^ + 14440640 x^ 
~ 6309980 x^ - 73169703 x^° + 282707420 - 61145056 

- 44575908x1* - 149321200x1^ - 446364 x^ - 2110212 x^. 

Appendix B.5. The order four and five operators for the staircase polygons. 

Let us give the order four and five operators -E5 and for Zq the staircase polygons 
generating functions [87] : 

_ 4 (3 - 140 x + 1295x2 - 1350x3) 
'"^^ x(l-x)(l-9x)(l-25x) 



(7 - 518 X + 6501 x2 - 8550 x^) 
^ x2 • (1-x) (l-9x)(l-25x) 
(1 - 196 X + 3963 x2 - 7200 x^) 
^ x3(l-x)(l-9x)(l-25x) 
1 - 57X+ 180x2 
x3 (1 - x) (1 - 9x) (1 - 25x)' 
p5 , in 1 - 70 x + 1176x2 -4032x3 

" x(l-4x)(l-16x)(l-36x) " 
25 - 2408 X + 51196 x2 - 211968 x^ 



2 



x2 (1 - 4 x) (1 - 16 x) (1 - 36 x) 
5 - 812 x + 23992x2 - 126720x3 ^ 
^ x3(l-4x)(l-16x)(l-36x) ' - 

1 - 516 x + 25956x2 - 193536x3 ^ 
^ x4 (1 - 4x) (1 - 16x) (1 - 36x) ' 
1 - 170 x + 2304x2 
x4(l-4x)(l-16x)(l-36x)' ^ ' ' 

Appendix C. Exponents of Fuchsian linear ODEs are generically algebraic 
numbers, not rational numbers 

Appendix C.l. "Lattice" Fuchsian ODEs. 

Keeping in mind some mainstream |127| conformal theory prejudice, the fact that 
critical exponents (for the ferromagnetic/ antiferromagnetic critical points) are rational 
numbers is too often taken for granted. However, here, we have a much stronger 
result: in all our previous calculations, the critical exponents all appear to be rational 
numbers for all the singularities of these ODEs (in our n-fold integrals of the Ising 
class, the exponents are even half-integers). Considering the polynomial coefficients in 
front of the successive derivatives in all our Fuchsian linear ODEs, one could wrongly 
imagine that this rational number exponent result is a straight consequence of the 
fact that these polynomial have, themselves, integer coefficients, this property being 
straightforwardly inherited from the enumerative combinatorics nature of the lattice 
problem. In that respect the simplest example of a Fuchsian linear differential operator 
certainly corresponds to the Gauss hypergeometric second-order differential operator 
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{D^ denotes the derivative d/dx) 

x-{\-x)-Dl + {c-{a + b+l)-x)-D,^ -ab, (C.l) 

which has the following indicial polynomials for the regular singular points x — 
0, 1, oo: 

r-(r-l + c), r-{r+ab/c), r-{r + l-c). (C.2) 

On these indicial polynomials one sees clearly that exponents being rational numbers 
is straightforwardly inherited from the rational character of the coefficients of (jC.l[) . 
More generally for Fuchsian linear ODEs of arbitrary order it can easily be shown, for a 
singular point which is a rational numbeiQ and occurring in the ODE with multiplicity 
one, that its exponents are necessarily rational numbers. 

However, it is important to underline that Fuchsian linear ODEs with integer 
coefficients do not have necessarily rational number exponents. Generically exponents 
of such integer coefficients Fuchsian linear ODE are algebraic numbers (algebraic over 
Q) not rational numbers. 

Appendix C.2. "Lattice" Fuchsian ODEs with algebraic numbers but not rational 
exponents. 

Let us consider, for instance, the order-four linear differential operator: 

{x^l)^x^ ■ A-{l-2x){l~x)-x- Dl (C.3) 

+ {l-6x + 6x^)-Dl+ 6-{l-2x)-D^ +12. 

This linear differential operator is Fuchsian. It has three regular singular points x = 0, 
X = 1 and X = oo, some of its exponents being algebraic numbers simply expressed 
in terms of the golden number as can be seen on the indicial polynomial corresponding 
respectively to a; = 0, 1, oo: 

r • (r — 1) • (r^ — r — 1), r • (r — 1) • (r^ — r — 1), 

(r-3) (r-2) (r + f)(r + 2), 

as well as on the solutions of this linear operator: 

(a; - 1)- V^+1/2 . ^1/2+1/2 VE . - 3 + 6 x) , 

^-1/2. V^+1/2 _ 1)1/2+1/2 . (_3 + 6 X + %/5) , (C.4) 

(3 a; — 1) • a:, 1 — 6 x^. 

It is interesting to calculate the p-curvature of this Fuchsian operator with a 
rational wronskian l/{{x — 1)^ • x'*), but with non-rational critical exponents (which 
cannot, therefore, be globally nilpotent). For an infinite number of primes the p- 
curvature is not nilpotent (a fortiori zero). However, for a subset of prime numbers (for 
instance, p = 11, 19, 29, 31, 41, 59, • • • , 281, 311, • • •) one finds a zero p-curvature. 
A heuristic interpretation is that, for some primes such that the golden number is 
"like" a rational numbe Jttl . an expression, like the first two ones in ljC.4[) , can be seen 
as an algebraic one. Of course this is not true for almost all primes. 

J And we have many examples 5 , in s or to or x, of rational number singular point, namely 
w = ±1/4, w = ±1/2, w = ±1,° X = 1/16, X = 1, X = 1/4, x = 1/9, x = 1/25, x = 1/8, . . . 
ft Namely the primes p such that — r — I factorize in Fp. For instance — r — 1 factorizes into 
(r + 7) (r + 3) mod. 11, (r + 4) (r + 14) mod. 19, ... , (r + 58) (r + 252) mod. 311. 
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Similarly, the order four linear differential operator 

(1 - x)2 • L*^ + 4 • (1 - 2 x) {1 - x) ■ x ■ (C.5) 

+ (l-Ux+lAx^) ■ Dl - 2- (l-2x) 

is also Fuchsian with solutions that can be expressed in terms of hypergeometric 
functions 3F2 and has quadratic number exponents as can be seen on the indicial 
polynomials corresponding respectively to x — 0, 1, coincide with l|C.4[) and for 
a; = 00 with (r — 1)^. Again calculating the p curvature of this Fuchsian operator 
with the same rational wronskian l/{{x — 1)^ • x"*), with algebraic but not rational 
numbers critical exponents (thus excluding global nilpotence) one finds, for a subset of 
the prime numbers, a nilpotent characteristic polynomial, namely (with a minimal 
polynomial T^). Generically the characteristic polynomial rules out the nilpotence, 
since it reads: 

+ {p^2)x3p +x2p ' ^ 

Another simple example is the order three linear differential operator 

■ {I - x) ■ Dl + X ■ {2 - 3x) ■ Dl + {l + 2x)-D^ -1, (C.7) 

which has exponents that can simply be expressed in terms of third root of unity and 
golden number, as can be seen on the indicial polynomial corresponding respectively 
to a; = 0, 1, 00: 

r-(r^-r + l), r-(r-l)^, r-(r^-3r + l). 

Again calculating the p-curvature of this Fuchsian operator with a rational 
wronskian l/((x — 1) • x^) but non rational critical exponents (thus excluding global 
nilpotence), one finds for a subset of the prime numbers (• • • , 109, 163, 181, 199, • • •) 
a characteristic polynomial with a minimal polynomial T^, and for a smaller set of 
primes (• • • , 73, 271, • • •) a characteristic polynomial with a minimal polynomial 

The miscellaneous examples we have displayed actually correspond to the generic 
situation of Fuchsian linear ODEs with integer coefficients (the proper framework we 
expect, at first sight, for lattice statistical mechanics quantities satisfying a linear 
ODE). However they are not globally nilpotent and thus are not "derived from 
geometry": genericalljl^ a Fuchsian linear ODE does not have solutions that can 
be expressed as n-fold integrals of algebraic integrands. 

Appendix D. Linear differential equation for Xd^'(0- 



decomposition 



The linear differential operator for v-) is of order eight, and has the direct sum 



£W ^ , with: = a + ^M-"^' 



dt ^ t 



(4) _ n3 , (5^' + 6t-l) {3t^ + 6t'-2t-l) 

~ ' +(i + f)(t-i)i-^* + ii + t)tHt-if 



2{l + t) {t-l)t^' 



(D.l) 



^ In a mathematical perspective. In contrast physics seems to favour the DFG framework for the 
minimal order ODEs . . . 
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where L^**^ is an order-four linear differential operator with apparent singularities 
~ lOt + 1 — 0, that will not be displayed here. Introducing the order-one operator: 

dt {t+lf 
one get rid of these apparent singularities (desingularization) , and obtains an order- 

-(4: 

^4 



five Fuchsian linear differential operator Gi ■ L^^^ which, after simple conjugaisons, 



can be simply written as: 

{t - 1)3 {t + 1)2 . Gi • if = {t- 1)3 {t + 1)2 ■ 

+ (-9 + llt^ + 26 1) {t - if {t + l)t^ ■ Df (D.3) 
+ (31 + 172 + 126 - 140 i + 19) it-l)t^ ■ 
+ 2 (11 + 107 1'^ + 179 - 271 + 74t - 4) < • D'^ 
+ {2t* + 4:3t^ + 327 f - 199t + 19)t ■ Dt +3 (t+1)^ 

The linear differential operator of order three, L^^^ is actually equivalent to 
the symmetric square Syvn?{LE) of the second order operator corresponding to the 
complete elliptic integral E{x^/'^) (see ([SSl) ). This order-three linear differential 
operator L^^^ is therefore globally nilpotent. Actually, we have calculated its p- 
curvature of L^^'^ and found that the corresponding characteristic polynomial (or 
minimal polynomial) reads . 

The order- four linear differential operator is also globally nilpotent: we have 
calculated the p-curvature and found that the corresponding characteristic polynomial 
(or minimal polynomial) reads T^. For the moment we have not been able to write 
one of his four solutions as a 4^F^ hypergeometric function up to a pull-back (trying 
to generalize subsection (|6.2p ). 

(n) 

Appendix E. Revisiting the global nilpotence of for n = 3, 4, 6. 

Appendix E.l. Revisiting the global nilpotence of • 

The global nilpotence of can be understood from the factorisation of the 

corresponding linear differential operator which can be seen as the direct sum of an 
operator of order three and of : 

Dx®L^, where: ~ Z2 ■ Li, with: (E-1) 

1 d /(1-f 3x +42;2)(1 -f-2a;)(a;- 1)^ 



92 ■ ^2 = q2 ■ Dl + 2 ■ {1 + x) ■ qi ■ Dr, + i ■ qq, where: 



2 

2, 



q2 = x(l-x)(l-|-4a;)(l-|-2a;)(l-4a;)(l + 3a;-|-4a;2) (1-l-a;) -g 
Q2 = 3264 X* -I- 56 x^- 862 x^-f 3641x5 4- 1873 x'^ 

-^149x^-23x2 4-2, 
qi = 1253376 x^^ + 1330688 x" - 492800 x" + 1432064 x^^ 
3680288 x" + 1249562x1" - 1 192677 x^ - 1051887 x^ 

- 317269 x^ - 47698 x^ - 8120 x^ - 2801 x"* - 693 x^ 

- 50x2 -M5x-F2, 
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qo = 626688 a;i5 + 1237248x1"* + 237504 898720x12 

+ 3726900 + 3657589 + 1424484 + 315618 

+ 122103 a;^ + 24147 - 21786 - 14389 ~ 3444 x^ 

- 375a;2-9a; + 2. 

The order- two linear differential operator Z2 can be seen to be homomorph to Qe, 
defined in (PT|) . corresponding to the complete elliptic integral E{4:x): 

Z2-Wi ^ W2-Qe, (E.2) 

where Wi and W2 are two linear operators of order one. From the explicit expression 
of Wi one easily finds the following solution for Z2: 

{12 x^ + 7 x'^ + X - 2) ■ E{Ax) 
^ (1 + 3 X + 4x2) (1 - 4 x) (1 + 2 a;) (1 - 

{Mx'^ + 11 x^ + 6 x^ + 7 x + 2) ■ K{4:x) 



(l + 2x)(l + x) (l + 3.T + 4x2)(l-2:)x' 
Appendix E.2. Revisiting the global nilpotence of ^^^^ . 

The global nilpotence of can be understood from the factorisation of the 

corresponding linear differential operator which can be seen as the direct sum of a 
linear operator of order-three and of D^'. 

Dx (B L3, where: L3 = L2 ■ Mi, with: 

1 d . / 4a; - 1 
,(x-l)2. 

L2 = q2 ■ Dl + (x - 1) ■ qi ■ Dx + go, where: (E.3) 

(72 = (1 - 16 a;) (1 -4x) (1024a;3 + 28x2 -42a;+ 1) (1 -a;)2 -x, 



qi ^ 262144 x^ - 228608 x^ - 4496 x^ + 19420x3 



3088x2 ^ 125 a;- 2, 



90 



147456 x^ - 242624 x^ + 13376 x^ + 49864 x^ 



14530x2 + 961x + 2. 



The order-two linear differential operator 

Qe - + ^ + ^ . (E.4) 

X (1 — 16 x) X 

corresponding to the complete elliptic integral £'(4-y/x), is equivalent to the linear 
differential operator L2 

Qe - zi = si • L2, (E.5) 

where Zi and Si are two order-one linear differential operators. 

Appendix E.3. Revisiting the global nilpotence of <1'^^. 

The linear differential operator for is an order-five Fuchsian linear operator which 
is the direct sum of Dx (here x = ?i;2) and of an order- four operator which factorises 
as a product of two order-two operators 

L.ie.) = Dx®La, where: L4 = M2 ■ L2, (E.6) 
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where M2 is a pretty large order-two linear differential operator (with a rational 
wronskian) and 

L2 = Dl - 2-{l -4x)-^-D, ^ 2^, with: 

P2 = (l-4a;)2(l-2:)(l-9a;)(l-10a; + 29a;2) 

X (1722 - 3306 + 2973 x* - 1548 x^ + 403 - 46 a; + 2), 

Pi = 898884 2797104 + 49021 17 X*- 5573337 a;^ 
+ 3999969 x^ - 1764005 x^ + 477136 x'^ - 79113 x^ 
+ 7883a;^-441x-l-ll, (E.7) 

Po = 898884 - 2559756 + 3491100 x^ - 2205501 x^ + 556746 x^ 
+ 92091 x^ - 92841 x"^ + 23740 x^ - 3081 x^ + 226 x - 8. 

The square of the wronskian of L2 is a simple rational function. The Fuchsian linear 
operator L2 is such that the p-curvatures are zero for almost all primes, and therefore 
it has a basis of algebraic solutions. Note that the differential Galois group of L2 
is isomorphic to the group of quaternions (eight elements). Its algebraic solutions 
correspond to an algebraic curve of genus g = 5. The equation of that algebraic 
curve reads: 

21025 • (2 - 36 X + 218 x^ - 558 x^ + 553 x'^ - 106 x^ + 27 x^)'' 
- 58 -^2 • (1 - 4x)2 (1 - x)2 (1 - 9x)2 (1 - 10x + 29x^)2 • 
+ (l-4x)''(l-x)'^(l-9x)'^(l-10x + 29x2)'^-Z^ = 0, (E.8) 

with: 

P2 ^ 1053 x^2 + 46836 x" - 429262 x^° + 520760 x^ + 1315505 x^ 
- 3318300 x^ + 3056140 x^ - 1518520 x^ + 448000 x"^ - 80280 x^ 
+ 8552x^-496x+ 12. 

Again, these algebraic functions, roots of a genus five algebraic curve, can be 
expressed as linear combinations of complete elliptic integrals of the third kind with 
a "characteristic" (first argument of the complete elliptic integral of the third kind) 
associated with a genus three curve (see [Appendix FT| ). 

We have also calculated the p-curvature of the (quite large) order-two Fuchsian 
linear differential operator M2 (for primes < 400) and found that all these p- 
curvatures are nilpotent. One can actually prove that M2 is equivalent to the previous 
second order operator (|E.4p . associated with E{4y/x). 

Appendix F. Towards a geometrical interpretation of global nilpotence 

Appendix F.l. Towards an interpretation as periods of algebraic varieties: closed 
formula for <i>^'' 

in] 

The integrals can all be expressed as sums of complete elliptic integrals of the 
third kind n{y{w), w), where the characteristic^ y = y{w) corresponds to some 

^ The first argument in a complete elliptic integral of the third kind is called the characteristic. 
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algebraic curves: 

^D^ = J2Mw)■Tl{y^iw),w), Pn{y^,w) = 0, (F.l) 

i 

where Ai (w) are algebraic expressions w and P„ are simple polynomials of yi and 
w with integer coefficients. 

Appendix F.1.1. Towards an interpretation as periods of algebraic varieties: closed 
formula for Let us give an exact expression for the integral (see (ITS)) ') in 

terms of complete elliptic integrals of the third kin41. Let us introduce /i and /2 
solutions of: 

(l-.T-4a;2)2.(/2 -1) -2x2.(4x + l)2.(/i-l) 

-2(l + 2.x)-(4a;-l)-(/i + l) - 0, (F.2) 

(l + 3x + 4a;2) . (/2_ 1) + (1 + 2 .x) (1 - 4 x) • (/a + 1) 

- (l + 4a;)2.(/2-l) = 0, (F.3) 
which are rational curves that can be parametrized as follows: 

Jl — 7. ^7 J2 — 



2(8 -u2)' (u2-6m-8)2' u2_gy_,_lg 

Let us introduce the three involutive birational transformations: 



The non-involutive birational transformation I o J maps (|F.2[) onto (|F.3p . and, of 
course, Jo/ maps back (IF.3[) onto (|F.2p . These two rational curves are invariant 
under the (Hadamard) involution H: fi l//i and /2 1//2- Note that H and 
/ o J commute. 

Let us introduce the expression: 



Rif, 4.) . . (IfiZ^mi:!^, (F.4) 

TT V / 



one then easily finds that i?(/, 4x) is such that 
Rif, 4x) + i?(l//. Ax) ^ 




{4xf -l){4x -f) 



f 

^-R{f{x), 4x) = Q{x) ---Eiix) + P{x) ■ - ■ K{Ax) 



(F.5) 



dx ■ TT ' TT 

for some (3(a;) and P{x) that can be deduced from / — y{x). These identities are 
valid for any y{x) and do not require y{x) to be, for instance, a rational function of 
X. The occurence of several operators equivalent to Le (associated to the complete 
elliptic integral E) in the factorisation of the linear differential operators corresponding 
to the ^^'^ , can be seen as a consequence of that identity (IF.5P on the complete elliptic 
integral of the third kind for arbitrary characteristic y{x). 

tt In the following subsections the calculations are expressed in terms of a variable x that is equal to 
w for n odd and to w'^ for n even. 

X We thank M. Rybowicz for kindly providing to us other closed expressions in terms of complete 
elliptic integrals of the third kind. 
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Introducing a combination of four complete elliptic integrals of the third kind, 
with characteristics satisfying the rational curves (|F.2p . (jF.3[) : 



E = i?(/i, Ax) - Ax) + R{f2, Ax) - i?(l//2, Ax), (F.6) 

can then be simply written in terms E as: 

$(3) - i 1 {l + Ax){l-9x^-12x^) 

- s + 12^ {l^x~Ax^){l + 3x + Ax^) ^(4x) + C 



where: 



C= ' ^^^"^^ (F7) 

48 y/1 + 2 xy/T^^/1 + 3x + Ax^' 



Appendix F.1.2. Towards an interpretation as periods of algebraic varieties: closed 
formula for "I>^ 
rational curves: 



formula for '^'"^ . For , and in a similar way as for , one has now the two 



3-(l+4x)(/i + l)2+ (1 -4x)(/i-l)2 = 0, (F.8) 

27- (1 + 4x):e^(/2 + 1)'+ (1-4:e)(2-5x2)2(/2-1)2 = 0. (F.9) 

A simple parametrization reads: 

1 m2 + 3 _ l+u fl - u\ fu^-Au-%\'^ 

Ji — 1 J2 



A 1*2-3' l-it' \1 + u) Vu2 + 4u-9 

Let us introduce the three involutive birational transformations: 

{l-x^^_{l-Ax^ j,..^. 
J [JO - (i_4^2) _(i_:r2).j^' ^^^^ - j^' ^U.j - ^- 

the non-involutive birational transformation J o I maps the first curve (|F.8p onto the 
second curve (jF.9p . and of course its inverse I o J maps back (jF.9[) onto (|F.8p . Note 
that these two rational curves arc actually invariant under the (Hadamard) involution 
H: fl ^ l/fi and /2 1//2- The rational curve (|F.8|) has, of course, many other 
involutive birational automorphisms: 

Using the same form for E as in (|F.6[) . the solution <i>^'' can now be written simply 
as: 



8 a;'' +9 - 11x2 - 1 
+ 72. (2. +1) (.2 +3.-1) • + ^ ■ 



where: 
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Appendix F.1.3. Towards an interpretation as periods of algebraic varieties: closed 

,(5) $(6) ^(7) ^ ^(8) 



formula for '^'^d j '^^d '^'^'^ '^'^d ■ For higher values of n (n > 5) the algebraic 



curves Pn{yi, w) = 0, which occur in the closed formula (jF.ip for the <i>^'' are no 
longer genus zero curves but higher genus curves. For instance for the <i>^^ and 
these curves are genus three curves. In terms of g = Awy they read respectively: 
-Aq^ - 4(4^2- 3) -g^ + 16 (2 - 1) • g 

+ 8 + (4w^-3w^-iy + l) = 0. 

and: 

lOg^ -40g3 - 5(3u;^ ~2)-q^ +80 (Sw^ ~ l) ■ q 

+ 32 (u> - 1) (2 w + 1) (2 w - 1) (w + 1) = 0. 

For <i>]j and these curves are genus ten curves. 
For ^'^j , in terms of g = 4:wy, one gets: 

q^ -6q^ -6(-5 + 4w2)g4 + 16 (6 - 5) • 

+ 24 (6 w'' + 5 - 12 w2)g2 - 96 (w - 1) (w + 1) (3 - 1) • g 
+ 32 - 128 - 192 + 288 w* + 32 = 0. 
The vanishing conditions 6n = of the discriminants in q of these polynomials, 

(n) 

associated with , read respectively: 

^ (1 - 2 u;) (1 + 4 w) (1 - 2 w - 4 - 1)2 

X {1 +w -3u;^ ~4w^f = 0, 
^6 = (l + 4i«)(l + 2w)(l-2w)(l-4w)(l-10w2 + 29w*f = 0, 
= (1 +4w)(l-2iy -4w2)(l -2w -Su-^ + 8 w^) (1 - 'u;)^ 
{w^ -w^ -2w + l)^ -d^ = 
d7 = l + 3w-10w2-35w^ + 5w'' + 62w^ + 17w^-32w^-16u;®. 

With the exception of (1 +4 w) (1 — 4^) = in = 0, these polynomial conditions 
are, respectively, the singulariti 
w has been changed into {—w). 



are, respectively, the singularities of the linear ODEs for tf>^'', and where 



Appendix G. Atkin's modular curves and Weber's modular functions 

The classical modular curve |128j which corresponds to the duplication of the ratio of 
periods of the elliptic curves |129j j = j(r) j' = j(2r): 

+ 3 • 15^ • (16 - 4027j/ + 16 (G.l) 
- 12 • 30^ • {j + j') + 8 • 30^ = 0, 

of course symmetric by j <-^- j', is well-known to be a genus zero curve with the 
rational parameterization: 

(z + 16)3 (z + 256)3 ./2l^^„„. 

z \ z / 

This is the duplication formula of the Klein's absolute invariant. The involution 
z 2^'^ / z is the Atkin involution |89| . Recall that the rational variable z can 
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simply be expressed in term of the Dedekind eta function: 



.6 /^(2T)^24 



z = 2'-{^ ■ (G.3) 
The modular curve (jG.ip can also be parameterized as: 

^ - (l-/fc2)2fc4 ' J - (I_fc2)4fc2 ^ ^ 

where the occurrence of the Landen transformation becomes explicit, or : 

(1 - 16 + 16 w^f ., _ (256 w"* - 16 + 1) 



(1-I6w2)i(;8 ' ^ (l + 4u))2(l -4^t;)2^^;4■ 
Similarly the triplication formula of the Klein's absolute invariant can also be written 
rationally: 

(z + 27)(z + 3)3 ., (z + 27) (z + 243)3 , fi'^ 

3 = J3[z) = , J = = J3 — 

z z'' \ z 

where z — > 3^/z is, again, the Atkin involution. The elimination of the rational 
variable z yields the classical modular curve which corresponds to the triplication of 
the ratio of periods of the elliptic curves j — j{T) j' ~ j(3r): 

/ + f - f ■ j" + 2232 . f ■ {j 1069956 • [j f + j' f) 

+ 2587918086 + 36864000 • (j + /) • {f + 241433 jj' + j'^) 
+ 16777216000000 • (27 + 27 j'^ - 45946 j j') 
+ 1855425871872000000000- (j + f) = 0. (G.4) 

The genus zero classical modular curve (jG.4p is, of course, symmetric by j ^ j' . 
The rational variable z can simply be expressed in term of the Dedekind eta 
function pOlIm] : 

r-(^Y\ s = 6. (G.5) 

We will not write (though it is straightforward) the classical modular curve ([50)1 
which corresponds to j = j(t) — > j' — j(6r). This classical modular curve is again 
a genus zero curve, corresponding to a polynomial relation symmetric by j ^ j' 
with integer coefScients and it can simply be obtained from the elimination of z in 
the rational parameterization (j5ip with, again, an Atkin involution z — )■ 2^ • 'S'^ / z. 

Note that the rational functions occurring in the rational parametrization of all 
these genus zero classical modular curves are related together by some rational change 
of variables: 

. f z-{z + 8)\ . /zjz + 9)2x 
Jeiz) =^ J2( — = J3( -z — • (G.6) 



z + 9 

The integers N such that the modular curves, PN{j{T), j(A^r)) = 0, are genus 
zero is a highly selected set of integers corresponding to the Monstruous Moonshine 
phenomenon [1321 1133j . 
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